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1.  INTRODUCTION 


1.1  Program  Objectives 

The  timely  detection  and  short  range  forecast  of  regions  of  convec¬ 
tive  storms  which  are  hazardous  to  aircraft  operations  is  of  vital  impor¬ 
tance  to  the  safety  of  flight  and  the  efficient  management  of  the  nation's 
airspace.  Weather  radars  have  the  capacity  to  observe  the  fine  structure 
of  severe  weather  and  to  localize  regions  of  potential  hazard.  Crane 
(1976)  postulated  that  the  potentially  hazardous  regions  should  coincide 
with  the  boundaries  of  strong  updrafts  during  the  growing  stage  of  a 
convective  cell  and  that  these  regions  could  be  identified  by  radar- 
detected  local  reflectivity  maxima.  The  objective  of  this  research 
program  was  to  test  the  hypothesis  that  reflectivity  or  tangential  shear 
cells  observed  by  a  Doppler  weather  radar  would  be  associated  with  poten¬ 
tially  hazardous  turbulence. 

Specifically,  the  research  program  was  to: 

1)  develop  and  test  algorithms  for  the  fine-scale 
detection  of  severe  weather  areas, 

2)  develop  and  test  algorithms  for  the  tracking  and 
extrapolation  of  movement  of  the  detected  turbulent 
areas ,  and 

3)  develop  ways  to  utilize  the  detection  and  tracking 
algorithms  in  enroute  and  terminal  air  traffic 
control  operations. 

Data  from  both  conventional  (reflectivity  only)  and  Doppler  weather 
radars  were  to  be  used  in  the  research.  The  algorithms  were  to  be  tested 
using  observations  from  instrumented  aircraft  obtained  during  severe 
weather  penetration  flights. 

1.2  Summary  of  Results 

1.2.1  Significant  Findings 

The  following  results  were  obtained  from  an  analysis  of  data  from 


1 


15  penetration  flights  supplied  by  the  South  Dakota  School  of  Mines  and 
Technology  (SDSMT)  and  of  simultaneous  weather  radar  observations  pro¬ 
vided  by  the  National  Center  for  Atmospheric  Research  (NCAR);  the  data 
were  from  two  of  the  National  Hail  Research  Experiment  (NHRE)  case  studies, 
22  June  1976  and  22  July  1976. 

1)  Cell  detection  algorithms  were  developed  which  success¬ 
fully  operated  on  reflectivity  or  on  Doppler  velocity 
variation  (shear)  data. 

2)  Tracking  and  position  extrapolation  algorithms  were 
developed  which  successfully  combined  reflectivity 
and  shear  data  obtained  at  different  elevation  angles 
and  generated  a  number  of  volume  cell  attributes  which 
were  used  to  evaluate  cell  significance. 

3)  A  total  of  135  encounters  with  moderate  or  more  intense 
turbulence  were  recorded  during  the  15  penetration 
flights.  Only  89  (67  percent)  of  the  encounters 

were  within  storm  areas  (contoured  regions)  having 
reflectivity  values  greater  than  20  dBZ*;  110  (82 
percent)  were  within  the  7  dBZ  contour  areas. 

Noting  that  most  weather  radar  systems  have  insuf¬ 
ficient  sensitivity  to  make  observations  over 
large  areas  at  levels  less  than  20  dBZ,  the  maximum 
probability  of  detecting  a  turbulent  patch  employ¬ 
ing  such  a  radar  is  less  than  67  percent. 

4)  Significant  reflectivity  cells  were  highly  correlated 
with  aircraft  encounters  with  turbulence.  It  was 
found  that  82  percent  of  the  moderate  or  more  intense 
turbulence  encounters  were  within  6  km  (3.2  nm)  of  a 
significant  reflectivity  cell.  For  encounters  within 
6  km  (3.2  nm)  of  a  region  having  20  dBZ  or  higher 
reflectivities,  the  probability  of  detection  was  94 
percent. 


*Radar  reflectivity  factor  expressed  in  decibels  relative  to  a  value  of 
1  mm°/m^. 
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The  addition  of  Doppler  information  did  not 
appreciably  change  the  detection  probability.  For 
turbulence  encounters  within  6  km  (3.2  nm)  of  a  signi¬ 
ficant  tangential  shear  cell,  the  two  C-band  radars 
averaged  a  71  percent  detection  probability;  the  two 
X-band  radars  averaged  a  57  percent  detection  pro¬ 
bability.  For  encounters  within  6  km  (3.2  nm)  of  a 
region  having  20  dBZ  or  higher  reflectivities,  the 
probabilities  were  increased  to  81  and  65  percent 
respectively. 

The  intensity  of  turbulence  encountered  by  the  air¬ 
craft  showed  little  correlation  with  either  cell 
reflectivity  or  tangential  shear.  The  observed 
tangential  shear  values  depended  on  the  relative 
locations  of  the  radars  and  the  cells. 

The  false  alarm  rate  (probability  of  cell  not  being 
turbulent)  depended  critically  on  the  assumed  radius 
of  influence  of  the  cell  and  the  intensity  of  tur¬ 
bulence.  For  moderate  or  more  intense  turbulence  and  a 
6-km  (3.2  nm)  radius  of  influence,  the  false  alarm 
rate  was  8  percent  for  the  significant  reflectivity 
cells  observed  by  the  S-band  radar,  averaged  9  percent 
for  the  C-band  Doppler  radars,  and  6  percent  for  the 
X-band  Doppler  radars. 

The  intensity  of  encountered  turbulence  depended  on 
the  age  of  the  cell.  The  most  intense  turbulent 
patches  were  associated  with  cells  which  were  first 
detected  at  or  just  after  the  time  of  aircraft  pene¬ 
tration.  Turbulent  patches  encountered  outside 
regions  with  precipitation  targets  suitable  for  radar 
detection  usually  developed  detectable  targets  after 
the  aircraft  encounter.  Successful  short  range 
forecast  of  the  probable  positions  for  new  cell 
development  will  be  mandatory  for  a  viable  tactical 
hazard  avoidance  system. 


9)  The  spatial  organization  pattern  of  new  cell  develop¬ 
ment  observed  in  the  two  Colorado  hail  storms  was 
similar  to  the  pattern  consistently  observed  by  Crane 
and  Hardy  (1981)  during  a  three-year  measurement 
program  in  Kansas.  The  constancy  of  the  observed 
spacings  between  developing  cell  clusters  suggests 
that  short  term  forecast  procedures  can  be  developed. 

1.2.2  Algorithm  Development 

The  research  effort  described  in  this  report  was  the  Federal 
Aviation  Administration  (FAA)  sponsored  segment  of  a  larger  program  of 
automated  radar  cell  detection  and  tracking  algorithm  development  and 
application  undertaken  by  Environmental  Research  8  Technology,  Inc.  (ERT) 
for  the  U.S.  Air  Force  Geophysics  Laboratory  (AFGL),  for  the  Water  and 
Power  Resources  Service  (WPRS),  and  for  the  FAA.  The  detection  and 
tracking  algorithms  have  been  described  by  Crane  (1979a).  The  segment 
conducted  for  AFGL  included  the  basic  cell  detection  algorithm  develop¬ 
ment  for  use  with  either  reflectivity  or  single  Doppler  velocity  data  and 
the  application  of  the  algorithms  to  real-time  processing  using  a  mini¬ 
computer  (Crane,  1979b).  The  segment  performed  for  WPRS  included  the 
development  of  tracking  algorithms  (Gustafson,  1980)  and  the  application 
of  these  algorithms  to  a  climatological  analysis  of  precipitation  in 
western  Kansas  for  use  in  evaluating  weather  modification  activities 
(Crane  and  Hardy,  1981). 

The  segment  performed  for  the  FAA  was  conducted  under  contract  to 
WPRS.  This  report  constitutes  the  final  report  for  the  1978-1981  time 
period  covered  by  the  contract.  Under  FAA  sponsorship,  algorithms  were 
developed  to  employ  both  tangential  shear  and  reflectivity  data  in  the 
tracking  process.  The  algorithms  were  designed  to  accept  data  from  more 
than  one  radar  thus  solving  the  multiple  radar  networking  problem  for 
cell  tracking  and  hazard  detection. 

The  radial  velocity  data  from  a  Doppler  radar  were  processed  to 
emphasize  local  velocity  perturbations.  In  the  analysis  performed  on 
the  NHRE  data,  the  magnitude  of  the  tangential  shear  was  employed  to 
measure  the  velocity  perturbation  at  the  edge  of  an  updraft  region  and 
elsewhere.  The  cell  detection  algorithm  operated  on  a  scalar  data 
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field  and  any  scalar  measure  of  velocity  perturbation  could  have  been 
used;  measures  such  as  radial  shear,  magnitude  of  the  vector  shear,  or 
Doppler  velocity  spread  (second  moment)  were  possible.  The  magnitude  of 
the  tangential  shear  was  used  in  this  study  because  (1)  it  is  optimum 
for  the  observation  of  an  azimuthally  isotropic  mesocyclone,  (2)  it  was 
employed  in  the  Joint  Doppler  Operational  Project  (JDOP,  1979)  for  severe 
weather  detection,  and  (3)  it  was  readily  calculated  for  each  of  the 
Doppler  radars.  Although  second  moment  data  are  often  recommended  for 
use  in  hazard  detection,  they  were  not  employed  in  this  study  because 
they  were  not  stored  on  the  data  tapes  for  all  the  Doppler  radars. 

The  cell  tracking  algorithm  generated  a  number  of  volume  cell  attri¬ 
butes  that  could  be  used  to  assay  cell  significance.  For  application  to 
hazard  detection,  a  significant  cell  was  a  cell  presumed  to  have  a 
turbulent  patch  within  a  reasonable  radius  of  influence  of  the  cell. 
Initially,  Crane  (1976)  viewed  the  radar  cells  as  isolated  building 
blocks  in  the  development  of  convective  storms.  He  expected  that  each 
cell  would  be  associated  with  an  updraft  region  and  that  the  observed 
reflectivity  maxima  would  be  horizontally  displaced  from  regions  of 
maximum  turbulence  occurring  at  the  interface  between  an  updraft  and  a 
downdraft  by  distances  the  order  of  half  a  cell  diameter.  Experience 
gained  from  processing  the  Kansas  data  for  WPRS  revealed  a  far  more 
complex  structure  to  the  convective  process.  Crane  and  Hardy  (1981) 
reported  at  least  four  different  cell  types,  only  one  of  which  should  be 
associated  with  significant  levels  of  turbulence. 

The  important  results  from  the  Kansas  data  analysis  were  (1)  that 
clusters  of  cells,  not  the  individual  cells  themselves,  were  the  consis¬ 
tent  indicators  of  the  important  updraft  regions  which  are  potentially 
hazardous,  and  (2)  that  significant  cells  were  cells  having  a  high  degree 
of  apparent  vertical  development.  The  concept  of  significance  applied  to 
the  processing  of  Kansas  data  included  clustering  and,  for  isolated  cells, 
high  reflectivity  and  moderate  vertical  development  or  significant  vertical 
development  in  the  absence  of  high  reflectivity.  The  resultant  entity 
was  called  a  significant  cell  and  cluster  (SC). 

The  cluster  detection  concept  did  not  work  when  initially  applied 
to  the  NHRE  data  because  of  the  inherently  high  false  tangential  shear 
cell  detection  rate  produced  by  Doppler  velocity  observations  at  low 
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signal-to-noise  ratios.  Unfortunately,  the  solution  to  the  cluster 
analysis  problem  was  worked  out  after  the  final  computer  tracking  runs 
were  completed  and  cell  clusters  could  only  be  used  in  the  analysis  of 
data  from  the  22  July  case  study  day. 

Vertical  development  and  persistence  of  either  reflectivity  cells, 
combined  reflectivity  and  tangential  shear  cells,  or  tangential  shear 
cells  were  used  in  this  study  as  an  indicator  of  significance.  As 
applied,  a  cell  had  to  be  detected  on  at  least  three  elevation  angles 
(span  a  height  interval  of  3  km  (10  kft)  or  more)  in  a  volume  scan 
sequence  (90  seconds  or  more  depending  on  the  radar  and  raster  scan 
limits)  and  had  to  exist  for  at  least  two  volume  scans.  It  is  noted 
that  JDOP  also  employed  vertical  development  and  persistence  criteria  in 
the  detection  of  mesocyclones .  They  employed  a  3-km  (10  kft)  minimum 
height  interval  restriction  but  a  longer  age  requirement,  5  to  15  minutes. 

Manual  analysis  of  the  spatial  locations  of  the  significant  cells 
detected  using  the  NHRE  case  study  data  revealed  that  a  significant  cell 
generally  occurred  as  a  member  of  a  cluster.  However,  the  number  of 
significant  cells  exceed  the  number  of  clusters  because  some  clusters 
contained  more  than  one  significant  cell.  The  net  result  was  that  the 
false  alarm  rates  reported  in  this  study  are  higher  than  they  would  have 
been  if  the  analysis  could  have  been  performed  using  clusters  to  detect 
turbulent  patches. 

1.2.3  Analysis 

The  objective  of  the  research  conducted  under  the  sponsorship  of 
the  FAA  was  to  evaluate  the  cell  approach  to  thunderstorm  hazard  detec¬ 
tion.  Data  from  five  separate  radars,  the  then  conventional  NCAR  CP-2 
S-band  radar,  the  two  NCAR  C-band  Doppler  radars,  CP-3  and  CP-4,  and  the 
two  NOAA  Wave  Propagation  Laboratory  X-band  Doppler  radars,  NOAA-C  and 
NOAA-D,  were  used  independently  to  generate  significant  cell  tracks 
whose  locations  were,  in  turn,  compared  with  simultaneous  aircraft 
encounters  with  turbulence  (see  Section  3.1  for  a  discussion  of  the 
radar  systems).  The  aircraft  encounters  were  used  as  the  turbulent- 
truth  indicators  of  hazard  (see  Section  4  for  a  discussion  of  aircraft 
observations) . 

Detection  was  claimed  where  the  cell  centroid  location  (at  ground 
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level)  was  within  a  specified  radius  of  influence  of  the  encountered 
patch  of  turbulence.  If  the  segment  of  aircraft  track  enclosed  within  a 
|  circle  centered  at  the  cell  centroid  position  and  having  the  specified 

radius  of  influence  was  not  turbulent  at  a  specified  or  higher  intensity, 
a  false  alarm  was  recorded  for  that  intensity.  The  total  counts  of 
turbulent  patches,  detected  patches,  radar  cells  within  the  radius  of 
I  influence  of  the  aircraft  track,  and  radar  cells  associated  with  one  or 

more  turbulent  patches  were  used  to  calculate  the  probabilities  of 
detection  of  a  turbulent  patch  or  the  false  alarm  rate  for  a  radar  cell 
observation. 

f  The  observed  detection  probabilities  for  each  of  the  Doppler  radars 

and  for  the  use  of  three  different  reflectivity  thresholds  for  the  con¬ 
ventional  S-band  radar  are  listed  in  Table  1  for  moderate  or  more  intense 
turbulence  and  a  6-km  (3.2  nm)  radius  of  influence.  Table  1  also  lists 
I  the  false  alarm  rate  and  a  critical  success  index  (see  Section  5.3) 

defined  as  the  number  of  successes  (cells  associated  with  turbulence) 
divided  by  the  total  number  of  detected  cells  within  the  radius  of  influ¬ 
ence  of  the  aircraft  track  plus  the  number  of  turbulent  patches  not 
i  associated  with  a  radar  cell  (failures).  Similar  statistics  for  encoun¬ 

ters  with  severe  or  extreme  turbulence  levels  are  listed  in  Table  2. 

These  tables  also  show  the  changes  in  the  probability  of  detection  and  the 
critical  success  index  that  occur  by  restricting  the  area  included  in  the 
t  analysis  of  the  detection  of  turbulence  to  (1)  all  segments  of  the  pene¬ 

tration  flights,  (2)  those  segments  within  a  region  enclosed  by  a  4-km 
(2.2  nm)  buffer  zone  outside  the  7  dBZ  contour  boundaries,  and  (3)  those 
segments  within  4  km  (2.2  nm)  of  the  20  dBZ  contour  boundaries. 

*  Variations  in  the  probability  of  detection,  false  alarm  rate,  and 

critical  success  index  with  the  radius  of  influence  presumed  for  each 
cell  are  displayed  in  Figures  1  through  3  respectively.  The  plotted 
data  are  for  turbulence  encounters  within  the  20  dBZ  contour  regions 
plus  buffer  zones  surrounding  the  contour  which  are  2  km  (1.1  nm)  smaller 
than  the  radius  of  influence.  The  data  are  displayed  by  frequency  band; 
the  results  were  combined  for  the  two  C-band  Doppler  radars  and  for  the  two 
X-band  Doppler  radars.  These  results  are  the  best  that  can  be  expected 
^  for  the  cell  approach  to  hazard  detection  because  the  analysis  has  been 

restricted  to  the  areas  of  the  radar  displays  having  precipitation  tracers 
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Figure  2  Probability  of  a  false  turbulent  patch  detection  (FAR);  T-28 
aircraft  at  6  km  height  and  90  m/s  airspeed 
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CRITICAL  SUCCESS 
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observable  by  radar.  They  show  that  the  radius  of  influence  of  the 
cells  was  in  excess  of  4  km  (2.2  nm).  By  a  6-km  (3.2  nm)  radius,  the 
rate  of  increase  in  detection  probability  and  of  decrease  in  false  alarm 
rate  had  slowed  perceptibly.  The  use  of  too  large  a  radius  of  influence 
is  not  recommended  because  of  the  increase  in  false  alarm  rate  and 
attendent  decrease  in  the  critical  success  index  expected  for  aircraft 
flights  through  the  weaker,  widespread  areas  of  a  storm  surrounding  the 
intense  updraft  regions.  Such  regions  were  not  adequately  probed  in  this 
study  because  the  primary  reason  for  the  aircraft  penetration  flights 
was  the  study  of  hail  development  in  the  most  active  regions  of  the  storm. 

The  critical  success  index  (CSI)  was  used  as  a  convenient  way  to 
combine  the  two  parameters  of  interest,  probability  of  detection  and 
false  alarm  rate.  High  CSI  values  can  only  be  achieved  when  the  detec¬ 
tion  probability  is  high  and  the  false  alarm  rate  is  low.  Referring  to 
Table  2,  it  is  noted  that  although  the  probability  of  the  detection  of 
severe  turbulence  was  high,  as  much  as  100  percent  in  some  instances, 
the  CSI  values  were  uniformly  low  due  to  the  high  numbers  of  false 
alarms.  For  aviation  safety,  the  probability  of  detection  is  the  most 
important  factor.  A  viable  hazard  detection  system  must  also  have  a  low 
false  alarm  rate  to  be  acceptable  to  the  aviation  industry.  A  combined 
index  such  as  the  CSI  was  used  in  this  study  to  provide  a  single  measure 
of  both  factors. 

The  CSI  values  for  reflectivity  data  alone  were  not  significantly 
different  from  the  CSI  values  obtained  for  the  C-band  Doppler  data  using 
the  magnitude  of  the  tangential  shear  to  measure  velocity  perturbations 
(Figure  3).  The  CSI  values  declined  significantly  when  only  high  reflec¬ 
tivity  (>40  or  50  dBZ)  cells  were  used  for  the  detection  of  turbulence. 
This  result  is  consistent  with  the  earlier  analysis  of  Burnham  and  Lee 
(1969)  but  differs  from  the  results  presented  by  Barclay  (1968,  1974) 
which  suggested  that  severe  turbulence  occurs  only  within  5.4  nm  (10  km) 
of  reflectivities  50  dBZ  or  higher.  The  results  from  the  two  NHRE  case 
studies  show  that  more  than  28  percent  of  the  encounters  with  moderate 
turbulence  occurred  outside  the  10  km  (5.4  nm)  buffer  zone  surrounding  a 
50  dBZ  contour.  Two  encounters  with  severe  or  more  intense  turbulence 
occurred  outside  the  10  km  (5.4  nm)  buffer  zone  surrounding  a  50  dBZ 
contour  (15  percent  of  the  encounters  with  severe  or  more  intense  turbu¬ 
lence)  .  Burnham  and  Lee  observed  that  lower  reflectivities  could  be 


associated  with  severe  turbulence  and  recommended  that  storms  with 
maximum  reflectivities  in  excess  of  40  dBZ  should  be  avoided  and  that  pene¬ 
trations  should  not  be  made  within  20  run  (37  km)  of  the  core  of  such  a  storm. 

Lee  (1977)  recognized  the  lack  of  correlation  between  the  reflec¬ 
tivity  values  along  an  aircraft  track  and  the  levels  of  turbulence 
encountered  by  the  aircraft.  He  investigated  the  use  of  several  measures 
of  velocity  perturbations  for  the  detection  of  turbulence  and  found  that 
Doppler  spread  was  the  most  promising.  He  reported  a  95  percent  proba¬ 
bility  of  detection  of  moderate  or  more  intense  turbulence  for  velocity 
perturbations  within  2  km  (1.1  nm)  of  the  turbulent  patches  but  did  not 
present  results  for  higher  turbulent  intensities  or  other  measures  of 
velocity  perturbations.  The  C-band  data  presented  in  Figure  1  show  a 
similarly  high  probability  of  detection  for  tangential  shear  at  a  radius 
of  influence  of  6  km  (3.2  nm) .  At  2  km  (1.1  nm) ,  the  probability  of 
detection  of  severe  turbulence  was  only  50  percent  but,  if  a  threshold  of 
severe  turbulence  was  selected,  the  probability  of  detection  was  100  per¬ 
cent.  In  the  latter  case,  the  false  alarm  rate  was  also  high  resulting  in 
a  critical  success  index  of  less  than  10  percent. 

The  X-band  radar  observations  produced  lower  CSI  values  than  did  the 
C-band  observations.  The  X-band  data  were  affected  by  attenuation  and 
the  signal -to-noise  ratios  for  observations  near  the  aircraft  track  were 
lower  than  for  the  other  radars  if,  in  fact,  a  signal  could  be  detected. 

The  net  results  were  lower  probabilities  of  detection  fsee  Figure  1)  and 
lower  CSI  values  than  for  the  C-band  observations.  The  data  presented  in 
Figures  1-3  included  encounters  both  within  the  4-km  (2.2  nm)  buffer  zone 
and  within  the  20  dBZ  contour  region  as  measured  at  S-band.  Therefore, 
the  loss  of  detectability  due  to  the  effects  of  attenuation  has  not  been 
corrected  in  the  X-band  (or  the  C-band)  data. 

Lee  (1977)  employed  Doppler  spread  (spectrum  breadth)  data  in  his 
analysis  of  thunderstorm  hazard  detection.  Spectrum  breadth  was  chosen 
because  of  the  expectation  that  the  turbulent  velocity  fluctuations  will 
be  isotropic  on  a  scale  commensurate  with  the  radar  range,  beamwidth 
product  (cross  range  dimension)  at  the  measurement  location  while  the 
larger  scale  shear  would  be  anisotropic.  He  deprecated  the  use  of 
radial  (Doppler)  velocity  shear  because  of  the  expectation  that  aniso¬ 
tropy  would  cause  the  magnitude  of  the  shear  to  vary  with  viewing  direc¬ 
tion  and  significantly  reduce  the  correlation  between  the  intensity  of  the 
turbulence  and  the  magnitude  of  the  shear. 
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Lee  reported  poor  correlations  between  turbulent  intensities  and 
the  magnitudes  of  reflectivity,  reflectivity  gradient,  first  or  second 
spatial  derivatives  of  the  radial  velocity  (shear  and  gradient  of  shear), 
or  Doppler  spread.  He  reported  high  probabilities  of  detection  when  he 
used  occurrences  of  Doppler  spread  above  a  preset  threshold  to  indicate 
moderate  or  more  intense  turbulence.  The  plots  that  he  presented 
(Figures  4.1  to  4.3  of  Lee,  1977)  also  showed  relatively  good  agreement 
between  the  magnitude  of  the  vector  shear  and  the  turbulent  intensity  but 
no  statistical  comparison  was  presented. 

In  this  study  we  found  that  the  occurrence  of  a  velocity  perturba¬ 
tion  as  indicated  by  a  tangential  shear  cell  was  a  good  indicator  of 
moderate  turbulence.  We  also  found  that  the  magnitudes  of  the  tangential 
shear  values  were  not  correlated  either  with  simultaneous  observations 
from  the  other  Doppler  radars  or  with  the  intensity  of  turbulence,  a 
result  indicative  of  anisotropy  on  the  scale  of  the  separation  between 
adjacent  radar  beam  positions  at  a  given  range  during  an  azimuth  scan. 
Anisotropy  on  this  scale  should  also  affect  detection  using  Doppler 
spread,  a  result  suggestive  of  one  of  the  underlying  causes  of  the  lack 
of  correlation  between  the  intensity  of  turbulence  and  Doppler  spread 
noted  by  Lee. 

Neither  the  reflectivity  values  nor  tangential  shear  values  were 
good  indicators  of  the  intensity  of  turbulence.  The  presence  of  a 
reflectivity  or  tangential  shear  cell  (perturbation)  was  a  good  indica¬ 
tor  of  moderate  turbulence  but,  as  the  intensity  of  the  turbulence  increased 
above  the  moderate  level,  the  critical  success  index  decreased  due  to  the 
lack  of  correlation  between  the  magnitudes  of  the  intensity  of  the  tur¬ 
bulence  and  the  cell  parameters.  For  application  to  the  detection  of 
moderate  turbulence,  this  lack  of  correlation  is  not  important.  For  the 
detection  of  severe  intensity  turbulence  in  a  widespread  region  of  less 
severe  weather,  the  lack  of  correlation  may  be  important.  In  the  latter 
case,  additional  information  such  as  cell  age  or  growth  rate  may  have 
to  be  employed  to  provide  a  higher  critical  sucess  index  for  hazard  detection. 

The  penetration  flights  were  into  regions  of  rapid  development  and 
strong  updrafts,  regions  that  should  be  turbulent,  regions  normally  avoided 
by  commercial  aircraft.  Exploration  of  these  highly  turbulent  regions 
could  bias  the  results  toward  lower  false  alarm  rates  and  higher  critical 
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success  indicies  than  would  be  expected  for  flight  paths  through  the 
regions  near  the  storms  that  could  be  safely  traversed  by  general  or 
commercial  aviation  aircraft.  This  bias  is  unavoidable  when  data  from 
hail  studies  are  used  for  the  evaluation  of  turbulence  detection  schemes. 
In  an  attempt  to  minimize  the  bias,  the  entire  flight  path,  both  within 
the  storm  and  while  maneuvering  before  and  after  each  penetration,  was 
used  in  the  statistical  analysis.  A  full  evaluation  of  the  hazard  detec¬ 
tion  algorithms  will  require  flights  in  and  around  a  wide  range  of  storm 
types  through  regions  that  are  expected  to  be  free  of  turbulence  as  well 
as  regions  that  are  expected  to  be  turbulent. 

1.2.4  Air  Traffic  Control  Applications 

The  ultimate  goal  of  this  research  program  is  an  automated  thunder¬ 
storm  hazard  detection  and  short  range  forecast  system  for  use  by  enroute 
and  terminal  area  air  traffic  controllers.  The  cells  (and  clusters)  are 
entities  which  can  be  readily  identified  and  displayed  to  a  controller 
along  with  aircraft  position  data.  Track  velocities  are  maintained  for 
each  cell  and  cluster  and  may  be  used  to  extrapolate  cell  positions  to 
forecast  the  locations  of  potentially  hazardous  regions.  The  potentially 
hazardous  regions  may  be  displayed  as  circles  of  cell  radius  of  influence 
(6  km  (3.2  nm) ,  say)  about  the  forecast  cell  position.  The  radius  for 
potential  hazard  may  also  be  expanded  to  allow  for  forecast  uncertainties. 

The  forecast  position  of  the  potentially  hazardous  region  must  be 
available  for  display  to  the  controller.  The  controller  is  interested  in 
the  probable  locations  of  turbulence  along  the  flight  path  some  5  to  20 
minutes  ahead  of  an  aircraft  encounter,  not  in  turbulence  that  has  been 
experienced  by  the  aircraft.  A  short  range  forecast  is  mandatory  because 
of  the  delays  inherent  in  collecting  the  radar  data,  processing  the 
data,  transmitting  the  data  to  a  center  and  finally  displaying  the  data 
to  a  controller. 

It  is  not  sufficient  to  display  only  potentially  hazardous  regions 
within  a  20  dBZ  contour  or  forecast  by  a  translation  of  the  hazardous 
regions  along  the  cell  trajectories.  Cell  lifetimes  are  relatively 
short.  Crane  and  Hardy  (1981)  found  that  the  average  cell  lifetime  was 
less  than  13  minutes  and  the  lifetimes  of  the  significant  cells  was 
roughly  a  factor  of  two  longer.  A  significant  fraction  of  the  observed 
cells  will  not  be  present  after  20  minutes  from  initial  detection  or, 
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using  a  minimum  of  two  volume  scan  sequences  to  estimate  a  cell  velocity, 
10  minutes  after  display  with  an  observed  velocity  for  extrapolation. 

More  importantly,  the  intensity  of  the  turbulence  tends  to  decrease  with 
age  and  a  forecast  of  probable  locations  for  new  cell  development  is 
required  to  provide  hazard  warning  at  or  before  the  time  a  cell  (or 
turbulent  patch)  is  detectable  by  radar. 

Moderate-to-severe  turbulence  was  observed  in  clear  regions  (reflec¬ 
tivities  less  than  7  dBZ)  within  15  km  (8.1  nm)  of  the  7  dBZ  contoured 
region  boundaries  on  one  of  the  days  (22  June  1976).  Severe  turbulence 
patches  existed  out  to  more  than  15  km  (8.1  nm)  from  the  closest  40  dBZ 
boundary,  to  19  km  (10.3  nm)  from  the  closest  50  dBZ  contour  boundary, 
and  to  20  km  (10.8  nm)  from  the  closest  multicell  storm  core  with  a  peak 
reflectivity  above  40  dBZ.  The  severe  turbulence  in  clear  regions, 
regions  outside  the  surveillance  capability  of  Doppler  weather  radars, 
were  outside  the  buffer  region  employed  by  the  aviation  turbulence 
advisory  service  operated  in  Australia  in  1974  (Barclay,  1974)  but  within 
the  region  to  be  avoided  as  specified  by  FAA  Advisory  Circular  00-24A 
(FAA,  1978)  and  recommended  by  Burnham  and  Lee  (1969).  Doppler  or  conven 
tional  weather  radar  observations  cannot  be  used  directly  to  reduce  the 
size  of  the  region  to  be  avoided;  but  forecasts  based  on  cell  and  cluster 
location  and  movement  information  should  be  useful  in  significantly  reduc 
ing  the  size  of  the  area  to  be  avoided. 

Analyses  of  an  extensive  set  of  significant  cell  and  cluster  (SC) 
data  from  Kansas  and  of  limited  data  sets  from  Colorado,  Oklahoma,  and 
Virginia  have  revealed  the  existence  of  a  preferred  spacing  between  SC's. 
Forecasts  based  on  the  preferred  spacing  concept  together  with  informa¬ 
tion  on  the  organization  and  motion  of  the  regions  of  severe  convection 
(detected  SC's)  should  be  developed  to  reduce  the  size  of  the  region  to 
be  avoided,  an  important  goal  of  hazard  detection  and  forecast  systems 
especially  for  terminal  area  operations.  What  is  evident  is  that  new 
cell  development  location  forecasts  are  required  to  solve  the  hazard 
detection  problem. 

Ultimately,  a  completely  objective  forecast  system  will  be  required 
to  generate  automatically  timely  hazard  warnings  for  display  to  a  con¬ 
troller.  Practically,  this  goal  will  not  be  reached  in  the  foreseeable 
future.  Our  current  knowledge  of  the  mesoscale  dynamics  of  storm  systems 
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is  imperfect  and  complete  automation  at  this  time  will  produce  more 
problems  than  solutions.  The  cell  detection  and  tracking  algorithms 
developed  during  this  research  program  and  the  allied  programs  for  WPRS 
and  AFGL  provide  the  essential  information  required  for  display  to  a 
meteorologist  for  manual  analysis  and  forecast.  Automated  forecast  pro¬ 
ducts  should  be  displayed  to  the  meteorologist  for  interpretation  prior 
to  display  to  a  controller.  These  products  should  also  be  recorded  for 
post  mission  analysis  and  forecast  verification.  With  experience,  the 
gradual  automation  of  the  entire  system  should  be  possible. 

1.3  Organization  of  Report 

This  report  summarizes  the  results  of  the  1978-1981  research  program 
conducted  for  the  Federal  Aviation  Administration  under  contract  to  the 
Water  and  Power  Resources  Service  of  the  Department  of  Interior.  Section 
2  reviews  the  background  of  the  thunderstorm  hazard  detection  problem 
and  the  application  of  the  cell  association  hypothesis  for  the  solution 
to  the  problem.  Section  3  presents  the  analysis  of  the  radar  observa¬ 
tions,  Section  4  describes  the  aircraft  penetration  data  and  Section  5 
provides  a  detailed  statistical  analysis  of  the  validity  of  the  cell, 
turbulence  association  hypothesis.  Recommendations  for  further  research 
are  presented  in  Section  6. 


18 


2 .  BACKGROUND 


2.1  Thunderstorm  Hazards 

Thunderstorms  have  long  been  known  to  be  potentially  hazardous  to 
aircraft  operations.  Several  types  of  hazards  exist.  Hail  encounters 
can  present  major  difficulties.  As  an  example,  the  Southern  Airways  DC-9 
crash  on  April  4,  1977  was  attributed  to  the  ingestion  of  massive  amounts 
of  water  and  hail  during  a  penetration  of  an  area  of  severe  weather 
(AWST,  1978).  Turbulence  and  low  level  wind  shear  encounters  also 
present  safety  hazards.  Brunstein  (1971)  reported  58  air  carrier 
accidents  during  the  1965-1969  time  period  attributed  to  turbulence  pro¬ 
duced  by  convective  activity.  Fujita  (1978)  documented  a  number  of  recent 
air  carrier  accidents  attributed  to  low  level  wind  shear  induced  by  cold 
outflow  from  convective  showers .  The  microburst  outflows  affect  aircraft 
while  landing  or  taking  off  and  therefore  are  important  to  the  terminal 
area  controller. 

Occurrences  of  hail  or  heavy  rain  are  readily  associated  with  local 
increases  in  reflectivity,  i.e.  reflectivity  cells,  observable  by  weather 
radar.  The  association  between  turbulence  in  convective  storms  and 
radar  observables  is  more  tenuous,  especially  if  only  conventional  radar 
data  (reflectivity)  are  available.  An  association  between  a  microburst 
detectable  by  direct  observations  only  near  the  ground  and  either 
reflectivity  or  Doppler  radar  observations  at  heights  well  above  the 
ground  is  even  more  tenuous. 

2.2  The  Cell  Approach  to  Hazard  Detection 

Crane  (1976)  argued  that  the  cloud-to-smal 1  mesoscale  structure  of 
the  reflectivity  field  (1-10  km  scale  size)  should  provide  adequate  clues 
for  the  detection  of  thunderstorm  hazards  due  to  hail,  to  heavy  rain,  or 
to  convective  turbulence.  The  argument  was  based  on  the  physics  of  cloud 
scale  motions  within  the  active  regions  of  convection.  During  the  early 
stage  of  the  development  of  a  convective  cell,  an  updraft  driven  by  the 
increased  buoyancy  produced  by  the  release  of  the  latent  heat  of  conden¬ 
sation  transports  moisture  aloft  which,  when  condensed  and  converted  to 
ice,  is  detectable  by  radar  as  a  local  reflectivity  maxima.  The  stages 


of  development  of  an  isolated  cell  are  depicted  schematically  in  Figures 
4a-d.  The  initial  stage,  presented  in  Figure  4a,  depicts  the  updraft  air 
flow  and  the  initial  echo  region  which  is  detected  as  a  three-dimensional 
reflectivity  (volume)  cell.  Turbulence  is  expected  to  occur  in  the 
localized  region  of  higher  shear  at  the  edge  of  the  updraft.  In  this 
picture,  precipitation  tracers  are  not  available  in  the  region  of  turbu¬ 
lence  for  observation  by  a  Doppler  radar  and  hazard  detection  is  possible 
only  by  the  presumed  proximity  between  the  updraft  edge  and  the  volume 
(reflectivity)  cell.  The  turbulence  patch  is  depicted  as  a  cross-hatched 
region  along  the  aircraft  track.  In  this  and  following  figures,  the 
cross-hatched  region  represents  moderate  turbulence;  a  solid  black  region 
represents  severe  levels  of  turbulence  (see  Figure  5  and  Section  4). 

As  the  cell  develops,  more  precipitation  particles  are  deposited  by 
the  updraft  which  settle  and  produce  a  localized  downdraft  by  evaporative 
cooling  and  atmospheric  drag.  The  downdraft  concentrates  the  return  flow 
from  the  updraft  and  produces  a  higher  shear  and  more  turbulence  at  the 
updraft,  downdraft  boundary  (see  Figure  4b).  In  this  case,  precipitation 
tracers  are  available  within  the  turbulent  patch  and  the  velocity  pertur¬ 
bations  produced  by  the  updraft,  downdraft  combination  are  detectable  by 
Doppler  radar.  The  cell  detection  algorithms  developed  by  Crane  (1979a) 
define  either  reflectivity  cells  which  enclose  the  area  within  3  dB  of  a 
local  reflectivity  maxima  or  tangential  shear  cells  which  enclose  the 
region  within  1.5  m/s/km  (0.9  kts/1000')  of  a  local  maxima  in  the  magnitude 
of  the  tangential  shear.  Using  these  criteria,  a  displacement  must  be 
expected  between  the  centroid  locations  of  the  reflectivity  cell  and  the 
tangential  shear  cell  at  the  updraft,  downdraft  interface  associated  with 
the  cell.  This  displacement  is  depicted  in  Figure  4b  by  an  offset  between 
the  volume  cell  and  the  turbulent  patch.  The  volume  cell  is  a  three- 
dimensional  assemblage  of  associated  reflectivity  (and/or  tangential 
shear)  cells  detected  at  successive  azimuth  scans  in  a  volume  scan  (tilt) 
sequence . 

The  downdraft  may  accelerate  and  produce  a  microburst  when  it  and 
the  accompanying  rain  shaft  strike  the  ground  (Figure  4c).  In  this  case, 
a  volume  cell,  a  middle  level  turbulent  patch,  and  a  microburst  (ground 
level)  are  all  in  close  proximity  and  lie  within  the  radius  of  influence 
of  the  volume  cell.  Successive  observations  of  newly  developing  volume 
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Figure  4a  Stages  in  the  development  of  an  isolated  convective  cell 


21 


SINGLE  CELL- SECOND  STAGE 
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SINGLE  CELL-RAIN  OR  GROUND 


Figure  4c  Stages  in  the  development  of  an  isolated  convective  cell 


SINGLE  CELL- DECAYING  STAG 


Figure  4d  Stages  in  the  development  of  an  isolated  convective  cell 


cells  associated  with  isolated  updraft  regions  show  initial  cell  develop¬ 
ment  aloft,  followed  by  a  rise  in  cell  height,  then  a  sinking  of  the 
cell  (Crane  and  Hardy,  1981).  The  sink  rate  approximates  the  expected 
fall  velocity  of  a  pulse  of  precipitation.  From  the  first  few  observa¬ 
tions  of  a  cell,  the  time  at  which  the  rain  pulse  strikes  the  ground  may 
be  calculated  and  used  to  locate  possible  microbursts.  Not  every  cell 
will  produce  a  microburst,  however,  because  mixing  may  occur  between  the 
downdraft  and  the  environmental  wind  both  above  and  near  the  surface; 
this  mixing  will  reduce  the  magnitude  and  velocity  of  the  outflow. 

Finally,  the  updraft  weakens,  the  remaining  precipitation  debris 
settles,  clumping  into  smaller  cells  and  producing  reflectivity  maxima 
(volume  cells)  in  a  horizontally  inhomogeneous  bright  band  by  aggregation 
and  melting.  In  this  case,  depicted  in  Figure  4d,  one  or  more  volume 
cells  may  be  detected  but  no  turbulence  is  expected.  The  volume  cells 
appear  in  a  cluster  which  is  not  associated  with  turbulence. 

This  physical  model  as  initially  proposed  by  Crane  predicted  a  high 
degree  of  association  between  reflectivity  cells,  hail  in  the  cores  of 
the  reflectivity  cells,  turbulence  at  middle  levels,  heavy  rain  on  the 
ground,  and  low  level  wind  shear  produced  by  microbursts.  It  was  the 
basis  for  the  turbulence  hazard  detection  hypothesis  that  initiated  this 
research  program.  Concomitant  research  for  WPRS  revealed  that  in  one  of 
the  Kansas  1978  aircraft  case  studies  the  mode  of  cell  development 
switched  from  isolated  cells  as  depicted  in  Figure  4  to  the  multi-cell 
structure  depicted  in  Figure  5  (see  Section  4.4,  Crane  and  Hardy,  1981). 
In  this  case,  the  updraft  continued  to  strengthen  and  cells  formed  in  a 
cluster.  At  the  time  of  cluster  development,  the  aircraft  crews  noted  a 
shift  to  more  vigorous  development  and  a  stronger  updraft.  The  updraft 
region  coincided  with  a  weaker  reflectivity  region  surrounded  by  four 
volume  cells.  The  higher  cell  depicted  schematically  in  the  overhang 
region  in  Figure  5  appeared  with  a  high  initial  reflectivity  value, 
greater  than  40  dBZ.  As  illustrated  in  Figure  5,  the  cluster  of  cells 
is  to  be  associated  with  stronger,  perhaps  extreme,  turbulence.  If  hail 
is  present,  the  volume  cells  in  the  cluster  will  coincide  with  the 
individual  hail  shafts  if  the  radar  system  has  sufficient  resolution  to 
separate  the  shafts. 

The  multi-volume-cell  cluster  surrounding  a  strong,  and  persistent 


MULTIPLE  CELL-  (CLUSTER) 


Figure  5  Volume  cells  in  a  cluster 
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updraft  can  be  identified  with  individual  weak  echo  regions  in  multi-cell 
storms  or  with  a  supercell  (Browning  and  Foote,  1976).  In  the  detailed 
Kansas  case  study,  the  cluster  was  an  elen  ‘nt  of  a  multi-cell  storm  with 
separations  between  isolated  significant  cells  and  the  observed  cluster 
of  the  order  of  10  km;  the  average  spacings  between  isolated  significant 
cells  and  clusters  (SC's)  during  the  most  active  period  of  the  entire 
storm  was  11.7  km.  The  model  to  be  evaluated  therefore  evolved  into  a 
consideration  of  significant  cells  and  clusters. 

2.3  The  Use  of  Doppler  Radar  Data 

Historically,  attempts  have  been  made  to  identify  thunderstorm  tur¬ 
bulence  hazards  on  the  basis  of  reflectivity  or  Doppler  data  field 
features  that  could  be  matched  on  a  point-by-point  basis  with  turbulence 
encounters  during  aircraft  penetrations.  In  these  analyses,  features, 
such  as  the  maximum  reflectivity  value  or  reflectivity  gradient  value 
encountered  by  an  aircraft  when  its  track  was  superimposed  upon  the 
reflectivity  field  measured  at  the  height  of  aircraft  penetration,  were 
used  for  comparison  with  the  spatially  coincident  aircraft  observations 
of  turbulence.  Lee  (1977)  produced  similar  point-by-point  comparisons 
with  features  of  the  radial  velocity  (Doppler  velocity)  field  and  with 
the  Doppler  spectrum  width  field  after  finding  little  correlation  between 
the  point-by-point  reflectivity  features  and  aircraft  turbulence.  He 
found  a  significant  number  of  coincidences  between  Doppler  spread  values 
in  excess  of  a  threshold  value  of  4  m/s  and  moderate  or  more  intense 
turbulence.  Using  occurrences  of  Doppler  spread  values  in  excess  of 
4  m/s  as  predictors  of  moderate  or  stronger  turbulence,  he  obtained  a  93 
percent  probability  of  detection.  If  he  used  Doppler  velocity  spread 
values  observed  within  2  km  of  the  aircraft  track,  he  increased  the 
probability  of  detection  to  95  percent. 

In  contrast  to  point-by-point  matching,  the  cell  detection  model 
anticipates  coincidence  between  the  cell  and  turbulence  only  within  a 
radius  of  influence  of  the  cell.  The  coincidence  is  expected  over  a 
wide  range  of  heights,  from  ground  level  (microbursts)  to  the  upper 
levels  of  the  troposphere;  the  detection  or  delineation  of  regions  of 
potential  hazard  therefore  apply  to  all  heights  within  the  troposphere. 

In  proposing  the  cell  detection  hypothesis.  Crane  (1976)  contended  that. 
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by  relaxing  the  coincidence  requirement,  conventional  weather  radar  data 
would  be  as  useful  as  Doppler  radar  data  for  the  detection  of  thunder¬ 
storm  turbulence.  He  further  anticipated  that  Doppler  radars  could  not 
be  used  directly  for  the  estimation  of  the  intensity  of  turbulence  because 
of  anisotropy  on  the  scale  of  the  radar  resolution  volume.  Lee  (1977) 
confirmed  this  expectation  by  observing  a  lack  of  correlation  between  the 
magnitudes  of  the  Doppler  spread  and  the  derived  gust  velocity  values. 

Crane  (Figure  A1 ,  1976)  summarized  turbulent  velocity  measurements 
made  by  Sinclair  (1974),  measurements  reproduced  in  Figure  6,  which  show 
that  the  vertical  velocity  variance  produced  by  thunderstorm  turbulence 
is  anisotropic  at  scales  larger  than  about  200  m  (departure  from  the 
k  ^  ^  line  at  scale  sizes  larger  than  200  m  is  indicative  of  anisotropy, 
see  Monin  and  Yaglom,  1967).  Doppler  spectral  broadening  is  produced  by 
turbulence  induced  velocity  fluctuations  on  scales  commensurate  with  the 
maximum  dimension  of  the  radar  resolution  volume.  For  radar  observations 
with  maximum  resolution  volume  dimensions  in  excess  of  500  m  (at  ranges 
in  excess  of  30  km  for  a  1°  beamwidth  radar)  the  Doppler  spread  is 
produced  by  anisotropic  turbulence.  When  the  turbulence  is  anisotropic, 
the  observed  Doppler  spread  values  no  longer  depend  on  the  intensity  of 
the  turbulence  in  the  isotropic  region  (the  region  of  each  spectrum  at 
scales  marked  "range  for  maximum  aircraft  response"  and  smaller  scales) 
which  can  be  specified  by  an  eddy  dissipation  rate,  e,  but  depend  on  the 
departure  from  isotropy  produced  by  the  three-dimensional  cloud  scale 
updraft,  downdraft  motion  and  larger  scale  vertical  and  horizontal  shear. 
The  result  is  an  observed  Doppler  velocity  variance  which  depends  on  the 
azimuthal  location  of  the  radar  relative  to  the  turbulent  patch  but  does 
not  depend  on  e . 

The  magnitude  of  the  velocity  fluctuations  (perturbations)  increase 
as  the  scale  size  increases  up  to  scales  of  the  order  of  5  km  (cloud 
scale);  as  shown  in  Figure  6,  for  larger  scales,  the  vertical  velocity 
(and  corresponding  horizontal  velocity)  variance  does  not  continue  to 
increase.  Observations  in  the  1  to  10  km  size  range  are  optimum  for  the 
observation  of  velocity  perturbations  induced  by  turbulence.  Gate-to-gate 
(radial)  shear  or  beam-to-beam  (tangential)  shear  may  be  readily  observed 
on  this  scale  producing  larger  perturbation  signatures  than  can  be 
expected  from  Doppler  spread  observations.  Tangential  shear  appeared  to 
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be  optimum  for  the  observation  of  the  cloud  scale  velocity  perturbations 
because  of  the  success  reported  by  Burgess  (1975)  and  by  Brown  and  Lemon 
(1976)  in  using  tangential  shear  for  the  detection  of  tornadoes  and  meso- 
cyclones.  They  found  significant  correlations  between  tangential  shear 
peaks  (mesocyclones  or  tornadic  vortex  signatures)  greater  than  5  m/s/km 
and  occurrences  of  tornadoes,  wind,  and  hail.  Using  mesocyclone  and 
tornado  vortex  signatures  to  detect  severe  weather  or  tornadoes,  the  JDOP 
report  (1979)  indicates  a  probability  of  detection  of  70  percent,  a  false 
alarm  rate  of  24  percent  and  a  critical  sucess  index  of  60  percent  for 
the  1977-1978  seasons. 

The  justification  for  the  use  of  tangential  shear  for  severe  weather 
analysis  is  based  on  the  idealized  Rankin  vortex  model  (Figure  7)  for  the 
flow  within  a  tornado  or  mesocyclone  (Donaldson,  1970).  The  idealized 
vortex  is  azimuthally  symmetric  and  observation  should  not  depend  on  the 
location  of  the  vortex  relative  to  the  radar.  The  strong  updrafts 
experienced  in  multi-volume-cell  clusters  however  show  a  tendency  to 
block  the  flow  at  middle  levels  of  the  troposphere  producing  flow  patterns 
similar  to  the  one  illustrated  in  Figure  8.  In  this  case,  the  tangential 
shear  is  anisotropic  and  the  observed  magnitude  will  depend  on  the  direc¬ 
tion  to  the  radar  relative  to  the  direction  of  the  middle  level  flow.  It 
is  noted  from  the  magnitude  of  the  tangential  shear  plots  that  more  than 
one  tangential  shear  cell  may  be  associated  with  a  vortex  or  blocking 
updraft . 

In  a  recent  review  of  the  use  of  single  Doppler  weather  radars 
(Wilson  et  al.,  1980),  the  use  of  radial  shear  was  demonstrated  for  the 
detection  of  shear  lines.  Either  radial  shear  or  tangential  shear  may  be 
employed  to  detect  localized  or  linear  velocity  perturbation  features. 
Radial  shear  observations  cannot  be  used  for  the  detection  of  idealized 
Rankin  vortices  but  idealized  azimuthally  symmetric  vortices  do  not  exist 
and  real  vorticies  should  produce  measurable  radial  shears.  Similarly, 
idealized  horizontal  shear  regions  along  constant  radius  arcs  cannot  be 
detected  using  tangential  shear  observations  but,  again,  they  do  not 
occur  in  nature.  The  use  of  either  one  or  the  other  or  the  magnitude  of 
the  vector  sum  of  the  two  shears  depends  on  which  measure  of  velocity 
field  perturbation  is  least  susceptible  to  measurement  uncertainties. 

Crane  (1977)  investigated  the  use  of  radial  and  tangential  shear  in  the 
analysis  of  weak  convective  activity  in  New  England  and  found  that 


LIFTING  WITH  SUPPRESSED  ENTRAINMENT 


Figure  8  Mid-level  flow  pattern  for  a  blocking  updraft  (Reynolds 
Number  of  40) 
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tangential  shear  perturbations  (cells)  were  less  affected  by  the  statis¬ 
tical  variations  in  velocity  measurements  than  were  radial  shear  cells 
|  and  recommended  the  use  of  tangential  shear. 
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3.  RADAR  OBSERVATIONS 


3,1  Significant  Cells  and  Clusters 

The  cell  detection  and  tracking  algorithms  developed  by  Crane 
(1979a)  provided  the  significant  cell  output  used  for  comparison  to 
aircraft  encounters  with  turbulence  in  the  evaluation  of  the  cell 
approach  to  hazard  detection.  In  brief,  the  radar  data  were  (1)  read 
from  magnetic  tapes,  (2)  calibrated,  (3)  averaged  to  roughly  a  1.5  km 
range  interval,  (4)  contoured  to  establish  the  regions  of  data  to  be 
processed  through  the  cell  detection  program,  and  (5)  processed  through 
cell  detection;  the  cells  were  then  tracked  (associated  from  one  azimuth 
scan  to  the  next  to  generate  volume  cells  for  each  raster  scan  sequence). 
Data  from  the  Doppler  radars  were  used  to  prepare  reflectivity  and 
tangential  shear  maps  for  each  azimuth  scan  and  both  data  fields  were 
processed  independently  by  the  cell  detection  program.  The  output  from 
cell  detection  consisted  of  an  interleaved  array  of  reflectivity  and 
tangential  shear  attributes.  The  output  from  tracking  the  reflectivity 
and  tangential  shear  cells  included  volume  cell  centroid  location,  track 
velocity,  and  the  average  values  of  reflectivity  and  tangential  shear. 

3.1.1  Conventional  Radar  Data 

The  radar  data  employed  for  analysis  were  acquired  in  two  different 
formats  and  processed  using  two  different  computer  systems.  Data  from 
the  conventional  CP-2  S-band  radar  were  acquired  in  the  packed  DADS  7- 
track  radar  data  format*  after  being  reformatted  to  9-track  tape.  Table 
3  contains  a  listing  of  the  characteristics  of  the  radars  which  have 
provided  data  for  processing  through  the  cell  detection  and  tracking 
programs.  The  CP-2  data  were  read  and  processed  by  an  early  version  of 
the  cell  detection  program  which  operated  on  an  IBM  370/158  computer. 

The  program  provided  contour  line  data  (vectors)  and  reflectivity  cell 
attributes.  The  contour  lines  were  plotted  as  shown  in  Figure  9  together 
with  the  cell  centroid  locations.  Cell  attributes  were  recorded  on 
magnetic  tape  and  stored  for  later  use.  The  attributes  included 
reflectivity,  area,  centroid  location,  time,  and  elevation  angle  (height) 

♦supplied  by  R.  Rinehart,  NCAR 


TABLE  3 


Figure  9 


EAST  (km) 


Computer  generated  contour  and  cell  location  plot.  Data 
are  for  0.6°  elevation  angle,  1640:50  MDT  on  22  July  1976. 
The  contours  enclose  reflectivity  values  >40  dBZ. 
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Figure  9  presents  observations  obtained  on  the  lowest  elevation 
angle  azimuth  scan  of  a  raster  scan  sequence  on  22  July  1976.  Contour 
data  were  prepared  at  two  thresholds,  17  and  40  dBZ.  Only  the  40  dBZ 
contours  are  displayed.  The  +  symbols  mark  the  centroid  locations  of 
reflectivity  cells  detected  at  the  same  elevation  angle.  The  Greek 
symbols  designate  the  persistent  clusters  (Section  3.2.1).  No  editing 
has  been  performed  on  the  data  and  both  ground  clutter  (ring  at  about 
10  km)  and  precipitation  targets  are  displayed.  The  track  for  aircraft 
penetration  Number  2  and  the  Greek  symbols  were  manually  superimposed  on 
the  figure.  Figures  10a  and  b  present  an  enlarged  view  of  the  area  around 
the  aircraft  track  at  the  height  and  time  of  penetration.  The  constant 
altitude  contour  plot  (Figure  10a)  was  provided  by  NCAR.  The  cell 
locations  and  intensities  displayed  in  Figure  10b  were  manually  obtained 
from  each  of  the  azimuth  scans  in  the  raster  scan  sequence  for  cell 
detections  within  1  km  in  height  of  the  6  km  height  constant  altitude 
plot. 

The  cell  detection  program  operated  on  all  the  input  data  corres¬ 
ponding  to  reflectivity  levels  greater  than  the  lowest  threshold  (17  dBZ). 
The  cells  enclosed  areas  with  reflectivity  values  within  3  dB  of  a  local 
maxima.  Only  areas  which  did  not  simultaneously  enclose  any  higher  level 
reflectivity  cells  were  used  to  define  a  cell.  Crane  and  Hardy  (1981) 

found  that  cells  defined  in  this  manner  were  relatively  small  with  an 

2 

average  area  of  5  km  .  The  cells  with  reflectivities  greater  than 
50  dBZ  displayed  in  Figure  10b  correspond  to  the  centers  of  the  50  dBZ 
contours  displayed  in  Figure  10a.  The  locations  of  the  other  cells  in 
Figure  10b  cannot  be  readily  compared  to  the  contours  in  10a  because 
of  the  differences  in  the  contouring  intervals  used  for  display  and  for 
cell  detection.  Of  interest  is  the  apparent  clustering  of  cells  within 
the  two  45  dBZ  contours  labeled  a  and  v.  The  significant  cell  and 
cluster  association  hypothesis  would  designate  the  region  of  aircraft 
track  between  1640  and  1641  flDT  as  having  a  high  probability  of  encounter¬ 
ing  severe  turbulence. 

3.1.2  Doppler  Radar  Data 

Doppler  data  from  the  NCAR  CP-3  and  CP-4  radars  and  the  N0AA  Wave 
Propagation  Laboratory  NOAA-C  and  NOAA-D  radars  were  supplied*  in  the 
♦supplied  by  F.I.  Harris,  NCAR 
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Detailed  contour  and  cell  locations  at  an  altitude  of  <>  km 
(msl)  for  1640  MDT,  22  July  1976.  Distances  are  measured 
relative  to  the  location  of  the  CP-2  radar  (Grover,  Colorado) 


expanded  FOFRAW  data  format  then  used  as  a  universal  format  for  NHRE  data 
storage  (see  Table  3  for  radar  characteristics).  These  data  were  read 
and  processed  on  the  CYBER-74  computer  at  the  Denver  Federal  Center  by  a 
modified  version  of  the  cell  detection  program  developed  for  WPRS  (see 
Gustafson,  1980  for  a  description  of  the  basic  WPRS  version  of  the  cell 
detection  program).  The  basic  program  was  revised  to  accept  the  FOFRAW 
data,  perform  the  required  calibration  analysis,  resolve  automatically 
velocity  ambiguities  (Crane,  1978),  calculate  and  average  tangential 
shear  estimates,  and  generate  tangential  shear  cells  by  operating  on  the 
magnitude  of  the  tangential  shear  field  with  the  same  cell  detection 
subroutine  as  used  for  reflectivity  analysis. 

Both  reflectivity  and  tangential  shear  cells  were  generated  for 
regions  having  reflectivity  values  in  excess  of  the  lowest  processing 
threshold  (17  dBZ) .  The  attributes  recorded  for  each  cell  type  were 
identical:  reflectivity,  magnitude  of  the  tangential  shear,  area, 
centroid  location,  time,  and  elevation  angle.  Contour  lines  were  gener¬ 
ated  at  two  reflectivity  levels,  17  and  40  dBZ  and  recorded  on  magnetic 
tape  for  future  display  but  none  was  plotted. 

The  cell  detection  algorithm  operated  on  quantized  scalar  data 
fields.  Cells  are  regions  enclosed  within  areas  corresponding  to  recorded 
observations  which  fall  three  quantization  steps  below  the  highest 
enclosed  value.  For  reflectivity  data  processing,  the  quantization  step 
was  1  dB  and  the  cells  enclosed  regions  within  3  dB  of  the  cell  maximum. 

For  tangential  shear,  the  data  were  quantized  in  1/2  m/s/km  shear 
intervals  and  cells  enclosed  regions  within  1.5  m/s/km  of  the  largest 
magnitude  shear.  Attributes  were  calculated  using  data  from  the  radar 
resolution  elements  included  within  a  cell  boundary.  The  average 
reflectivity  value  was  calculated  from  the  reflectivities  at  the  range, 
azimuth  locations  within  a  tangential  shear  cell;  conversely,  the  average 
magnitude  of  the  tangential  shear  was  calculated  from  the  tangential 
shear  values  enclosed  within  a  reflectivity  cell. 

3.1.3  Volume  Cell  Tracking 

A  statistical,  best  match  association  process  was  used  in  the 

automated,  computer  tracking  of  volume  cells  (Crane,  1979a).  A  Lagrangian 

tracking  procedure  was  employed  which  shifted  the  location  of  a  cell*  from 

‘volume  cell  refers  to  the  3-dimensional  entity  being  tracked;  a  cell  refers 
to  the  output  of  the  cell  detection  program  obtained  at  a  fixed  elevation 
angle  on  a  single  azimuth  scan. 
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its  position  at  the  time  of  detection  to  the  position  it  would  have  had 
at  the  reference  time  for  each  volume  scan.  The  volume  cells  were 
initially  given  the  average  velocity  of  all  the  cells  being  tracked  dur¬ 
ing  the  prior  volume  scan  and  the  velocity  of  an  individual  volume  cell 
was  adjusted  from  its  initial  value  as  tracking  proceeded.  A  starting 
average  track  velocity  was  estimated  from  radiosonde  data  for  each  of 
the  case  study  days  although  the  precise  starting  value  was  not  critical 
to  the  successful  operation  of  the  tracking  program. 

When  Doppler  radar  data  were  available,  the  tracking  program 
operated  first  on  the  reflectivity  cells  from  an  azimuth  scan,  then  on 
the  tangential  shear  cells  for  the  same  azimuth  scan.  Data  of  either 
type  were  handled  in  the  same  way  although  cell  associations  with  a 
volume  cell  were  tallied  by  type.  A  number  of  additional  attributes  were 
calculated  for  each  volume  scan  (tilt  sequence)  including  volume  cell 
translation  velocity,  height,  and  the  identities  of  the  contour  and 
cluster  in  which  it  was  imbedded. 

Clusters  were  defined  on  the  basis  of  volume  cell  similarity  and 
proximity.  The  best  match  association  procedure  tallied  all  the  cells 
(reflectivity  or  tangential  shear)  and  all  the  volume  cells  that  could 
be  matched  with  each  other  on  a  volume  scan.  The  best  cell,  volume  cell 
pairings  were  used  for  tracking  and  the  list  of  possible  matches  was 
used  to  identify  the  elements  of  a  cluster.  The  maximum  cell,  volume 
cell  separation  allowed  in  a  cluster  at  the  reference  time  was  4.5  km. 

A  cluster  with  a  number  of  volume  cells  could  span  a  larger  region  since 
a  cell  could  provide  the  bridge  between  two  volume  cells  and  a  volume 
cell  could  bridge  two  cells  spreading  the  domain  of  a  cluster  over  a 
number  of  volume  cells. 

Clusters  defined  in  this  manner  are  generally  small  in  horizonal 
extent.  The  maximum  cluster  size  observed  in  the  CP-2  (reflectivity 

2 

only)  data  from  the  22  July  1976  case  study  was  18  km  (area  of  82  km  ); 

the  average  maximum  dimension  of  a  best  fit  ellipsoid  spanning  a  cluster 

2 

was  5.8±.3  km  (average  area  13±1  km  ).  The  association  logic  was 
sensitive,  however,  to  false  cell  detections  which  were  common  in  the 
tangential  shear  data  at  low  signal-to-noise  ratios.  As  an  example, 
cluster  observations  at  1640  MDT  on  22  July  1976  are  presented  in  Figure 
11.  The  Doppler  radars  were  operated  with  a  five-minute  volume  (tilt) 
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scan  sequence.  The  data  from  the  conventional,  CP-2  radar  were  combined 
in  a  composite  volume  scan  spanning  the  same  five-minute  time  interval. 
Because  of  the  sequence  of  elevation  angles  employed  by  the  CP-2  radar 
during  a  five-minute  interval,  the  version  of  the  tracking  program  used 
for  analysis  reported  the  data  as  if  from  seven  volume  scans,  not  the 
expected  composite  volume  scan  nor  the  two  to  three  volume  scans  actually 
made  during  the  interval.  Unfortunately,  all  seven  reported  scans  were 
labeled  with  the  same  identification  number  making  separation  difficult. 

As  a  result,  all  the  data  are  plotted  with  an  additional  code  to  identify 
the  number  of  volume  cells  in  a  cluster.  The  Doppler  data  all  correspond 
to  single  volume  scans  and  the  existence  of  three  or  more  volume  cells 
in  a  cluster. 

The  clusters  labeled  a  and  v  in  Figure  11,  correspond  to  the  major 
cellular  impulses  H  and  G  depicted  by  Wade  and  Foote  (1978)  and  by 
Heymsfield  et  al .  (1980)  for  storm  III.  The  cellular  impulses  correspond 
to  the  clusters  automatically  separated  by  the  tracking  algorithm  for 
all  but  the  NOAA-C  radar.  The  low  reflectivity  region  between  cluster 
X  to  the  North  and  the  combined  a,v  region  of  higher  reflectivity  (see 
Figures  9  and  10a  and  also  Jameson  and  Heymsfield,  1980)  was  populated 
by  low  reflectivity  tangential  shear  cells  which  bridged  the  gap  between 
the  separate  clusters.  These  false  cells  are  in  evidence  in  Figure  12. 

In  this  figure,  only  the  significant  cells  were  plotted. 

For  this  study,  a  significant  cell  consisted  of  three  or  more  cells 
(reflectivity  or  tangential  shear)  associated  with  a  volume  cell  during 
the  five-minute  composite  volume  scan.  With  the  exception  of  cell  detec¬ 
tions  by  the  NOAA-C  radar,  the  cells  appear  as  a  number  of  tightly  grouped 
clumps  corresponding  to  the  clusters  in  Figure  11.  As  before,  multiple 
cells  are  depicted  at  each  significant  cell  location  due  to  the  increased 
number  of  apparent  volume  scans  in  a  composite  volume  scan  interval  for 
the  CP-2  radar.  The  wider  area  of  false  alarms  contributed  to  the  large 
area  cluster  depicted  in  Figure  11  for  the  NOAA-C  radar.  Since  the 
centroid  location  for  the  cluster  was  computed  using  reflectivity  weighted 
volume  cell  positions,  the  cluster  is  positioned  between  the  three  high 

reflectivity  clusters  corresponding  to  a,  v  and  X.  The  area  associated 

2 

with  this  cluster  was  139  km  ,  more  than  10  times  the  average  observed 
using  CP-2. 
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Fictitiously  large  and  improperly  positioned  clusters  were  also 
detected  by  each  of  the  Doppler  radars  on  22  June  1976  and  both  the 
NOAA-C  and  CP-3  radars  occasionally  produced  large  clusters  for  the  22 
July  case  study.  Because  of  the  large  clusters  (high  false  alarm  rates 
for  tangential  shear  cells)  and,  more  importantly,  because  the  computer 
tape  containing  the  cell  detection  output  for  the  CP-2  data  for  the 
aircraft  penetration  time  periods  on  22  June  1976  was  lost,  the  cluster 
data  were  only  partially  analyzed  for  the  evaluation  of  hazard  detection. 
Because  the  significant  cell  data  were  readily  available  and  reliably 
positioned,  the  principal  statistical  analyses  were  performed  using  only 
the  significant  cell  data.  The  increase  in  false  cell  detections  near 
the  aircraft  track  will  bias  the  statistics  towards  higher  false  alarm 
rates  than  would  be  expected  for  the  cluster  data  alone.  Because  the 
aircraft  penetration  flights  were  generally  through  the  higher  reflec¬ 
tivity  regions  of  the  storm,  the  effect  of  the  bias  on  the  evaluation  of 
the  utility  of  the  Doppler  radar  should  be  small.  The  cluster  and 
significant  cell  results  listed  in  Tables  1  and  2  are  not  significantly 
different. 

3.2  The  NHRE  Case  Studies 

Radar  data  were  supplied  by  NCAR  and  aircraft  data  were  supplied  by 
South  Dakota  School  of  Mines  and  Technology  (SDSMT)  for  two  of  the 
National  Hail  Research  Experiment  (NHRE)  case  study  days,  22  June  1976 
and  22  July  1976.  An  analysis  of  the  former  day  was  published  by  Harris 
et  al.  (1978);  analyses  of  the  latter  by  Wade  and  Foote  (1978),  Kelley 
et  al.  (1978),  Heymsfield  et  al.  (1978),  Foote  et  al.  (1978),  Jameson 
and  Heymsfield  (1980),  and  Heymsfield  et  al.  (1980). 

The  two  case  studies  correspond  to  the  occurrence  of  multicell  hail 
storms  optimumly  positioned  for  triple  Doppler  radar  analyses  of  the 
three-dimensional  wind  fields.  Four  Doppler  radars  were  employed.  Their 
characteristics  have  been  summarized  by  Harris  et  al.  (1978)  and  by  Wade 
and  Foote  (1978).  The  important  radar  characteristics  are  listed  in 
Table  3.  Their  locations  are  depicted  in  Figures  11  and  12.  Extensive 
manual  (computer  aided)  editing  was  required  to  prepare  the  data  for 
multiple  Doppler  radar  analysis.  Because  of  the  large  manpower  commit¬ 
ment  required  for  multiple  Doppler  analysis,  only  a  limited  number  of 


volume  scans  were  prepared  for  each  of  the  case  studies.  Although  the 
automated  processing  programs  employed  in  this  research  did  not  require 
editing,  the  data  to  be  processed  for  the  case  studies  were  prepared  in 
the  universal  format  by  NCAR,  and  consequently  only  these  data  to  be 
edited  were  available  for  further  analysis.  As  a  result,  only  6  to  11 
volume  scans  (depending  on  the  radar)  or  less  than  55  min  of  data  were 
available  for  the  22  June  case  study  and  2  to  4  volume  scans  (less  than 
20  minutes)  for  the  22  July  case  study.  The  scans  were  spaced  by  inter¬ 
vals  as  large  as  15  to  20  minutes  negating  their  utility  for  automated 
tracking. 

Continuous,  rapid  volume  scan  (1.5  minutes  per  scan)  data  were 
available  for  the  CP-2  radar.  Because  of  the  missing  computer  output 
tape,  only  data  from  22  July  could  be  automatically  processed  but, 
because  of  the  elevation  angle  sequence  employed  by  the  radar,  the 
resultant  volume  cells  could  not  be  tracked  automatically  from  one  volume 
scan  to  the  next.  The  tracking  program  (and,  for  CP-2  and  22  June,  the 
plotted  data)  did  provide  significant  volume  cell  data  which  could  be 
readily  compared  to  the  aircraft  observations.  Figure  12  provides  a 
sample  of  such  a  comparison. 

Triple-Doppler  generated  velocity  fields  were  available  for  two  of 
the  volume  scans  for  22  June  (Harris  et  al„,  1978)  and  for  four  volume 
scans  for  22  July  (Kelly  et  al . ,  1978).  A  detailed  cell  and  storm 
development  analysis  was  also  available  for  the  latter  storm  (Wade  and 
Foote,  1978;  Foote  et  al.,  1978).  The  latter  storm  is  considered  in  more 
detail  in  this  report  because  of  the  wealth  of  material  for  comparison. 

Superficially,  both  case  study  days  were  similar.  A  number  of  cells 
-  major  (peak)  cellular  impulses  (PCI)  -  were  observed  which  developed  on 
the  upwind  flank  of  the  storm,  moved  through  the  storm,  and  decayed  on 
the  downwind  side.  This  type  of  storm  behavior  was  observed  in  Kansas 
(Crane  and  Hardy,  1981)  and  has  been  documented  many  times  (Renick,  1971; 
Marwitz,  1972).  However,  the  22  July  storm  revealed  an  abrupt  change  in 
its  mode  of  development  part  way  through  the  period  of  analysis  (Foote, 
et  al.,  1978).  Cluster  data  for  the  22  July  storm  provided  an  additional 
means  for  observing  the  change.  The  lack  of  sufficient  data  on  a  scale 
larger  than  that  provided  by  the  limited  surveillance  volumes  imposed  by 
the  rapid  raster  scan  observations  make  a  detailed  analysis  impossible, 
however. 


'I 


3.2.1  Storm  Structure,  22  July  1976 

Wade  and  Foote  (1978)  described  the  morphological  development  of 
Storm  III  on  22  July  in  terms  of  a  succession  of  major  convective 
impulses  which  developed  over  the  1500-1725  MDT  time  period  and  moved 
toward  the  south.  They  observed  the  impulses  in  the  CP-2  reflectivity 
data  near  storm  top  (11  km)  and  tracked  the  peak  reflectivity  impulses 
using  5  dBZ  contour  maps  at  a  constant,  7  km,  altitude.  The  volume 
cells  and  clusters  of  volume  cells  that  formed  during  the  same  time 
period  can  also  be  used  to  observe  and  record  storm  development  and  to 
compare  storm  descriptors.  Differences  are  to  be  expected  between  the 
two  methods  for  depicting  storm  structure.  Persistent  clusters  should 
be  indicative  of  active  regions  of  convection  accompanied  by  a  sustained 
updraft  and  associated  transient  updraft,  downdraft  couplets  and  volume 
cells.  Each  reflectivity  peak  used  by  Wade  and  Foote  should  correspond 
to  a  single  volume  cell  in  a  cluster  which  moves  through  the  cluster  as 
both  the  volume  cells  and  cluster  evolve. 

CP-2  data  from  1539  to  1722  MDT  were  processed  through  the  cell 
detection  and  tracking  programs.  The  resultant  cluster  detections  were 
plotted  for  four  successive  25-minute  time  intervals  (identified  as  I  to 
IV)  in  Figures  13a-d.  The  persistent  clusters  apparent  in  the  data 
were  labeled  by  Greek  letters  a  to  t.  Where  appropriate,  the  peak 
reflectivity  convective  impulse  (PCI)  identity  provided  by  Wade  and 
Foote  (1978)  was  included  as  well.  Cluster  a  (PCI  D  at  1540  MDT)  was 
tracked  from  one  volume  scan  to  the  next  and  its  average  velocity  was 
used  to  adjust  for  cluster  motion.  The  average  velocity  for  cluster  a 
during  its  observation  period,  1539-1652  MDT,  was  12°/4 .4  m/s,  which  was 
almost  identical  to  the  average  storm  velocity,  10°/4.3  m/s,  reported  by 
Wade  and  Foote  for  storm  III  for  the  1500  to  1630  time  period.  This 
velocity  was  used  for  the  Lagrangian  adjustment  of  the  cluster  positions 
in  each  of  the  four  plots.  Figures  13a-d. 

The  data  for  tracking  period  I  (Figure  13a)  show  tightly  grouped 
cluster  locations  each  corresponding  to  a  separate  persistent  cluster. 
Transient  clusters  seen  only  on  one  or  two  volume  scans  were  not  identi¬ 
fied  as  persistent.  The  remnants  of  Storm  I  and  clusters  from  Storm  II 
(cluster  6  plus  transient  clusters  perhaps  representing  the  decaying 
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Figure  13b  Locations  of  clusters  relative  to  cluster  a  for  four  successive  25-minute  time 
periods  on  22  July  1976.  Geographical  location  relative  to  Grover,  Colorado 
(CP-2)  can  be  obtained  from  the  plotted  locations  of  the  (100,-10)  coordinate 
for  each  of  the  5,  5-minute  periods  in  a  time  interval.  The  plotted  symbols 
are  identified  in  Figure  11. 
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state  of  an  earlier  persistent  cluster)  are  also  displayed.  Two  of  the 
major  convective  impulses,  D  and  E,  passed  through  cluster  a  during  this 
time  interval.  The  volume  cells  associated  with  each  impulse  moved 
apparently  toward  the  north  and  east  through  the  cluster  finally  produc¬ 
ing  the  clusters  reported  at  locations  about  5  km  from  the  origin.  The 
cluster  pair  at  6  km  initially  appeared  at  that  position  and  was  transient 

Each  of  the  clusters  displayed  a  relatively  small  spread  of  its 
constituent  parts,  a  result  indicative  of  little  variation  in  the  track 
velocities  from  one  cluster  to  the  next.  An  apparent  increase  in  the 
spread  of  a  cluster  is  usually  indicative  of  the  demise  of  one  cluster 
and  the  formation  of  a  new  cluster.  Cluster  8*  and  n*1  (the  superscript 
designates  the  time  period)  were  adjacent  to  each  other  (compare  Figures 
13a  and  b)  but  the  transition  from  one  to  the  next  was  made  in  a  discrete 
step,  hence  the  different  labels.  The  x's  in  cluster  n  (two  volume  cells 
in  a  cluster)  appeared  first  separated  by  about  10  km  from  the  dots  (•)  of 
cluster  8,  then  the  area  within  cluster  n  began  to  fill  in.  Similarly, 
cluster  0  formed  as  a  discrete  identity  separated  by  about  8  km  from  0 
and  10  km  from  n. 

New  clusters  formed  as  discrete  entities  along  an  east-west  line 

just  to  the  north  of  cluster  a  is  shown  by  clusters  a1,  B* ,  n**,  011  and 
III 

tt  in  Figure  14.  These  clusters  were  forming  on  the  upwind  flank  of 
storms  II  and  III  along  a  broad,  60-km  line.  During  time  period  II  a  new 
cluster  formed  to  the  north  (downwind)  of  cluster  a.  This  cluster  was 
separated  from  a  by  a  region  of  low  reflectivity  and  its  development  may 
be  characterized  as  anomalous  rather  than  normal.  New  cluster  develop¬ 
ment  then  shifted  to  the  northwest,  southeast  line  through  y*1*,  v  , 

,111  IV  ,  IV  „  ,  ..  ,  .  . 

C  ,  a  and  x  .  Development  along  this  line  was  normal  with  new 

IV 

clusters  forming  on  the  upwind  flank.  Cluster  p  ,  however,  developed 
downwind  of  5*^  along  a  line  parallel  to  the  x^  line. 

The  abrupt  shift  in  the  mode  of  cluster  development,  from  along  an 
east-west  line  to  along  the  y*^,  x line  was  noted  by  Foote  et  al. 

(1978).  It  appears  to  be  in  response  to  a  larger  scale  modulation  of 
regions  of  preferred  cell  development.  Such  modulation  has  been  noted 
by  Crane  and  Hardy  (1981).  This  case  differs  from  the  Kansas  observations 
because  the  pattern  presented  in  Figure  14  was  translating  relative  to  the 
terrain  at  12°/4.4  m/s  while  the  patterns  observed  in  Kansas  were 
stationary  relative  to  the  terrain. 
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These  data  also  showed  similarities  to  the  Kansas  observations  in 
the  average  values  of  the  horizontal  spacings  between  clusters.  For  the 
1629-1723  time  period  (III  and  IV),  the  average  spacing  between  SC's 
(clusters  plus  isolated  significant  cells)  was  14.8+.5  km,  and  for  the 
earlier  time  period  (I  and  II),  the  average  spacing  between  SC's  was 
13.4±.8  km.  By  way  of  comparison,  the  average  value  for  the  36  storm 
elements  in  the  Kansas  HIPLEX  sample  was  14.3H.0  km,  a  value  essentially 
identical  to  the  values  for  the  22  July  case  study. 

The  water  flux  per  SC,  a  value  found  by  Crane  and  Hardy  (1981)  to 
vary  little  from  one  convective  storm  to  the  next  was  within  2o  (two 
standard  deviations)  of  the  average  value  observed  in  Kansas  (after  correc¬ 
tion  for  the  known  difference  between  the  Kansas  radar  and  rain  gauge 
estimates  of  water  flux).  The  corrected  value  for  Kansas  was  0.44  million 
metric  tons  per  hour  per  SC  (MT/h/SC);  the  value  for  time  periods  III  and 
IV  for  the  22  July  1976  storm  was  1.31.3  MT/h/SC.  Crane  and  Hardy  showed 
that  the  relative  values  of  water  flux  per  SC  could  be  used  as  an  indica¬ 
tor  of  overall  storm  efficiency.  Barring  radar  calibration  differences, 
the  22  July  storm  could  be  judged  as  being  more  efficient  than  the  Kansas 
HIPLEX  sample.  The  Kansas  data  were  for  the  entire  25  to  150  km  annular 
surveillance  region  of  the  Skywater  75  radar  but  the  22  July  data  were 
from  the  more  limited  surveillance  area  of  a  raster  scan  imposed  by  the 
limited  coverage  of  a  dual  or  triple  Doppler  radar  system.  Noting  that 
observations  were  made  of  just  the  most  active  cells,  the  data  from 
Colorado  should  be  biased  toward  higher  values. 

The  NCAR  analyses  of  the  22  July  case  study  included  both  cell 
morphology  and  flow  field  depictions  constructed  from  the  multiple  Doppler 
radar  data.  The  horizontal  flow  fields  at  the  inflow  and  mid-levels  of 
storm  III  at  the  time  of  aircraft  penetration  Number  2,  1640  MDT,  are 
illustrated  in  Figure  15  (adapted  from  Kelley  et  al . ,  1978).  The  per¬ 
sistent  clusters  are  identified  at  both  levels.  The  detailed  reflectivity 
field  at  mid-level  (6  km  mol)  is  shown  in  Figure  10,  the  detailed  cluster 
field  is  presented  in  Figure  11  and  the  significant  cell  locations  are 
depicted  in  Figure  12.  Several  features  of  the  flow  field  are  apparent 
in  Figure  15:  blocking  at  low  levels  by  a  downdraft  in  cluster  v, 
blocking  by  the  updraft  in  cluster  v  at  mid-level  (6  km  mol),  and 
acceleration  of  the  flow  between  a  and  v.  In  each  instance,  velocity 
perturbations  are  evident  which  should  be  detectable  by  a  single  Doppler 
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Figure  15  Horizontal  flow  field  at  the  inflow  and  mid-level  of  storm  III  for  the 

1637-1641  MDT  time  period,  22  July  1976  (adapted  from  Kelley  et  al . ,  1978) 


radar.  The  mid-level  convergence  region  to  the  northwest  of  cluster  6 
should  also  be  detectable  as  well  as  the  mid-level  perturbations  in 
clusters  a,  6,  and  A  (due  to  the  persistent  updrafts  in  the  clusters; 
clusters  a  and  6  are  more  than  an  hour  old  at  the  time  of  Figure  15). 

An  examination  of  Figures  12  and  15  reveals  that  the  velocity 
perturbations  are  detectable  using  tangential  shear  to  observe  the 
features.  Each  of  the  Doppler  radars  detected  the  mid-level  perturba¬ 
tions  in  cluster  A,  v  and  a.  The  CP-4  radar  apparently  could  readily 
sense  the  low  level  blocked  flow  upwind  of  cluster  a.  The  extended 
region  of  velocity  perturbations  (acceleration)  observed  by  the  radar 
produced  a  number  of  closely  spaced  tangential  shear  cells  which  combine 
to  produce  a  cluster.  The  clusters  observed  by  the  other  Doppler  radars 
only  depicted  the  accelerated  flow  evident  higher  in  the  storm  and  between 
a  and  v  (Figure  11).  Each  of  the  radars  detected  6  as  a  cluster. 

The  updraft  cores  were  identified  and  marked  by  heavy  arrows 
corresponding  to  the  magnitude  and  direction  of  the  environmental  wind  at 
3.5  km  height  (inflow)  with  the  magnitude  increased  by  the  ratio  of  the 
density  of  the  air  at  3.5  km  to  the  density  at  6.5  km  height  (mid-level). 
If  the  air  is  lifted  without  entrainment  (see  Figure  5)  the  air  in  the 
updraft  core  should  conserve  horizontal  momentum  while  being  lifted  from 
one  height  to  the  other.  Using  this  simple  model,  the  updraft  cores 
could  be  readily  detected  by  either  multiple  or  single  Doppler  radar 
systems.  The  core  locations  depicted  in  Figure  15  correspond  to  the 
locations  of  the  mid-level  wind  vectors  having  the  appropriate  magnitude 
and  direction.  The  analysis  presented  by  Kelley  et  al.  (1978)  and  by 
Heymsfield  et  al.  (1980)  concentrated  on  cluster  v  (their  PCI  I)  corres¬ 
ponding  to  the  weak  echo  region  produced  by  the  updraft  to  the  south  and 
east  of  the  symbol  v  in  Figure  15b.  They  showed  a  significant  updraft 
extending  from  the  updraft  core  location  on  Figure  15b  through  the  region 
of  accelerated  flow  between  a  and  v.  The  aircraft  data  showed  the 
updraft  core  (highest  value  along  the  aircraft  track)  to  lie  just  to  the 
west  (within  2  km)  of  the  location  indicated  on  Figure  15b.  The  recent 
reanalysis  of  the  flow  field  reported  by  Heymsfield  et  al.  (1980) 
indicates  that  the  updraft  core  (observed  from  the  magnitude  and  direc¬ 
tion  of  the  mid-level  flow)  should  be  positioned  1.6  km  to  the  west  of 


the  location  indicated  in  Figure  15b  producing  even  better  agreement 
with  the  aircraft  measurements. 

3.2.2  Intercomparison  Between  Radars 

The  radars  employed  for  the  22  July  case  study  were  operated  to 
provide  overlapping  surveillance  of  the  region  with  strong  convection 
enclosing  clusters  a,  v,  X,  6,  £  and  a.  As  illustrated  in  Figure  12, 
each  of  the  significant  cells  were  observed  by  more  than  one  radar.  A 
listing  of  the  reflectivity  and  tangential  shear  reported  for  each 
significant  volume  cell  by  radar  is  given  in  Appendix  A  for  the  volume 
scans  which  coincided  with  the  seven  aircraft  penetration  flights. 

Doppler  data  were  provided  by  NCAR  for  four  of  the  seven  volume  scans 
and  Doppler  data  were  available  from  all  four  radars  for  only  two  volume 
scans. 

Figure  16  intercompares  the  S-band  (CP-2)  and  C-band  radar  data  for 
each  of  the  simultaneously  observed  volume  cells.  The  reflectivity 
values  used  in  the  analysis  were  average  values  for  the  volume  cell  as 
observed  by  the  radar.  A  comparison  between  the  CP-2  and  CP-3  observa¬ 
tions  show  agreement,  and  one  can  conclude  that  there  is  no  calibration 
differences  between  the  radars.  Comparison  between  CP-2  and  CP-4  or 
between  CP-3  and  CP-4,  as  shown  in  Figure  17,  show  a  decided  calibration 
error  for  CP-4.  The  median  difference  between  the  CP-4  reflectivity 
observations  and  the  other  two  radars  was  10  dB. 

No  correlation  was  evident  between  the  magnitude  of  the  tangential 
shear  and  reflectivity  or  between  the  magnitudes  of  the  tangential  shear 
values  observed  by  the  two  Doppler  radars,  CP-3  and  CP-4.  The  former 
result  was  expected  on  the  basis  of  the  physical  models  for  cell  develop¬ 
ment  (Section  2)  and  on  the  basis  of  prior  aircraft  and  radar  observa¬ 
tions  of  convective  storms.  The  latter  result  is  doe  to  the  anisotropy 
of  the  velocity  perturbations  (turbulence)  and  to  the  differences  in  the 
observing  scale  defined  by  the  separation  between  the  positions  of  the 
centers  of  the  resolution  elements  used  in  the  calculation  of  tangential 
shear  (Section  2).  Referring  to  Appendix  A  and  Figure  15  for  observa¬ 
tions  during  penetration  2,  the  velocity  perturbations  in  clusters  a,  v 
and  S  were  observable  by  all  four  radars.  The  tangential  shear  observed 
by  CP-4  for  cluster  v  ranged  from  one-half  to  three  times  the  nearly 
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Figure  17  Simultaneously  observed  reflectivity  and  tangential  shear  values 
comparison  between  the  two  C-band  radars,  22  July  1976 
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constant  value  observed  from  CP-3.  The  increased  variability  evident  for 
CP-4  is  most  likely  due  to  the  differences  in  range  between  the  volume 
cells  and  the  radar;  CP-4  was  about  15  km  from  cluster  v  while  CP-3  was 
about  55  km  from  the  same  cluster. 

Comparisons  between  (1)  the  X-band  and  S-band  radar  (Figure  18),  (2) 
the  two  X-band  radars  (Figure  19),  or  between  an  X-band  and  a  C-band 
radar  (Figure  20)  yielded  similar  results.  Again,  no  consistency  was 
evident  between  the  tangential  shear  observations  of  the  two  radars.  An 
apparent  calibration  difference  was  evident  but,  because  of  fundamental 
limitations*  in  the  X-band  recording  equipment,  the  difference  could  be 
due  to  the  effective  saturation  of  the  reported  X-band  received  power 
observations  and  the  differences  in  ranges  between  the  radars  and  the 
volume  cells.  Little  significance  was  attributed  to  the  X-band  reflec¬ 
tivity  data;  it  was  used  primarily  to  set  processing  thresholds. 

The  proximity  between  detections  of  a  volume  cell  by  different  radars 
and  the  small  differences  in  observed  reflectivity  between  well  calibrated 
radars,  less  than  5  dB  rms,  indicate  that  the  cell  detection  and  tracking 
algorithms  are  identifying  the  same  features  in  a  storm.  The  simultan¬ 
eous  observations  may  be  used  to  intercalibrate  the  radars.  Cell  detec¬ 
tion  data  from  several  radars  may  be  processed  simultaneously  by  the 
tracking  program  to  provide  a  single,  best  estimate  output  for  the  net¬ 
work  of  radars.  Refinements  are  possible  to  remove  differences  in  the 
sizes  of  the  resolution  elements  (differences  in  the  ranges  from  the 
radars  to  a  volume  cell)  and  to  combine  the  radial  velocity  data  to 
detect  the  positions  of  updraft  cores  in  a  more  realistic  manner. 


‘private  communication  with  J.  Miller  of  NCAR 
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Figure  18  Simultaneously  observed  reflectivity  and  tangential  shear  values 
comparison  between  the  X-band  and  S-band  radars,  22  July  1976 
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Figure  19  Simultaneously  observed  reflectivity  and  tangential  shear  values 
comparison  between  the  2  K-band  radars,  22  July  1976 


4.  AIRCRAFT  PENETRATIONS 


Aircraft  penetration  flights  were  made  through  the  most  active 
regions  of  convection  within  the  bounds  of  multiple  Doppler  radar  cover¬ 
age  on  each  of  the  case  study  days.  The  aircraft  employed  for  the  pene¬ 
trations  was  the  armored  T-28  operated  by  the  South  Dakota  School  of 
Mines  and  Technology  (Sand  and  Schleuseuer,  1974).  The  aircraft  was 
instrumented  to  record  the  thermodynamic  properties  of  the  air,  to 
observe  liquid  water  content,  hail  mass,  ice  crystal  habit  and  vertical 
acceleration,  and  to  calculate  turbulence  levels  and  vertical  velocity. 

For  this  study,  only  the  vertical  velocity,  vertical  acceleration,  and 
turbulence  levels  were  of  interest.  Recorded  values  for  penetration  2 
on  22  July  1976  are  presented  in  Figure  21.  These  data  correspond  to 
the  radar  observations  presented  in  Figures  9,  10,  11,  12,  13c  and  15. 

The  aircraft  track  was  located  in  the  figures  displaying  the  radar 
data  on  the  basis  of  smoothed  position  data  prepared  by  NCAR.  Several 
position  determination  systems  were  available  such  as  an  M-33  transponder 
tracking  radar  colocated  with  the  CP-2  radar,  onboard  distance  measure¬ 
ment  equipment  (DME) ,  and  the  FAA  Air  Traffic  Control  transponder  network. 
The  M-33  transponder  data  were  used  in  this  analysis.  With  smoothing, 
the  estimated  position  error  is  0.5  km*. 

The  vertical  velocity  values  displayed  in  Figure  21  were  calculated 
from  either  the  aircraft  rate  of  climb  or  from  observed  pressure  altitude 
changes  and  corrected  for  the  flight  characteristics  of  the  aircraft. 

The  expected  measurement  uncertainty  is  less  than  5  m/s*.  The  vertical 
acceleration  values  were  instantaneous  values  measured  at  the  center  of 
gravity  of  the  aircraft  which  were  sampled  once  per  second.  The  turbu¬ 
lence  levels  were  calculated  from  observed  true  air  speed  fluctuations 
(Sand  et  al . ,  1976).  The  true  air  speed  was  obtained  from  the  indicated 
air  speed  which  was  sampled  twice  per  second.  The  fluctuations  were 
spectrally  analyzed  (16  point  FFT)  to  produce  a  velocity  variance  value 
in  eight  different  frequency  (spatial  scale)  bands.  The  variance  values 
were  converted  to  eddy  dissipation  rate  (e)  estimates  and  averaged  to 
provide  a  best  estimate  value  of  e  consistent  with  an  inertial  subrange 
spectrum  (Monin  and  Yaglom,  1975).  The  results  were  reported  as  e1//3  in 

*P.  Smith  of  SDSMT,  private  communication 
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Figure  21  Vertical  velocity,  vertical  acceleration  (departure  from  quiet  -  IG) ,  and  (cduy  dissipation 
rate)1/3  observed  during  penetration  2,  1639-1644  MDT,  22  July  1976;  T-28  aircraft  at  6  km 
height  and  90  m/s  airspeed 
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units  of  cm  s  Based  on  the  spectra  presented  by  Sand  et  al . ,  the 
rms  uncertainty  in  is  estimated  to  be  25  percent. 

4.1  Turbulence  Levels 

Two  systems  have  been  used  to  characterize  atmospheric  turbulence, 
one  based  on  eddy  dissipation  rate  measurements  (MacCready,  1964)  and 
the  other  employing  derived  gust  velocities  (Pratt  and  Walker,  1954; 

Lee,  1974).  The  latter  measure  of  turbulence  is  proportional  to  the 
peak  vertical  acceleration  change  encountered  in  a  short  time  interval 
(order  of  a  second).  The  derived  gust  velocity  is  the  strength  of  a 
standardized  isolated  gust  that  would  have  produced  the  observed  vertical 
acceleration  change.  The  proportionality  coefficient  depends  upon  the 
parameters  of  the  aircraft  that  relate  to  the  calculation  of  lift.  The 
derived  gust  velocity  is  used  to  provide  an  estimate  of  the  level  of 
turbulence  that  does  not  depend  on  the  characteristics  of  the  observing 
aircraft.  It  may  be  used  to  estimate  the  response  (vertical  accelera¬ 
tion  changes)  of  other  aircraft  to  the  same  idealized  gust. 

Eddy  dissipation  rate  is  a  measure  of  the  turbulent  fluctuations 
of  the  air  that  is  independent  of  the  response  of  a  particular  aircraft 
to  the  turbulence.  MacCready  (1964)  showed  that  the  rms  value  of  vertical 
acceleration  (gust  loading)  was  linearly  related  to  e^3  and  suggested 
that  e^3  be  used  to  characterize  turbulence  because  it  does  not  depend 
n  the  characteristics  of  an  aircraft  and  it  can  be  used  to  determine  the 
rms  gust  loading  for  the  aircraft. 

Both  e^3  and  the  peak  acceleration  deviations  from  quiet  condi¬ 
tions  (1  g)  in  a  short  time  interval  were  recorded  during  the  T-28 

1/3 

penetration  flights.  These  data  were  used  to  compare  the  largest  e 
values  within  a  30-second  time  interval  (3  km  of  flight)  and  the  largest 
instantaneous  acceleration  departure  from  quiet  conditions  within  the 
same  time  interval.  A  one-sided  measurement  of  acceleration  was  employed 
rather  than  the  largest  instantaneous  change  used  in  the  estimation  of 
a  derived  gust  velocity  because  the  acceleration  observations  were  made 
twice  per  second  but  displayed  only  once  per  second  and  large  exersions 
were  noted  as  occasionally  occurring  between  the  recorded  samples.  The 
results  of  comparing  the  two  observations,  peak  acceleration  and  e^3, 
are  presented  in  Figure  22.  As  predicted  by  MacCready,  a  straight  line 
approximation  provides  a  good  fit  to  the  observations. 
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The  eddy  dissipation  rate  is  a  quantitative  measure  of  the  intensity 
of  the  turbulence  in  the  atmosphere.  The  more  familiar  intensity  descrip¬ 
tors,  light,  moderate,  severe,  and  extreme,  are  not  quantitative  but  are 
used  to  represent  the  subjective  response  of  a  pilot  to  an  encounter  with 
turbulence.  A  pilot  responds  to  rapid  changes  in  acceleration;  the  rms 
gust  loading  or  the  peak-to-peak  vertical  acceleration  changes  within 
the  time  interval  of  about  a  second  may  be  used  as  nearly  direct  measures 
of  the  reaction  of  a  pilot  to  the  turbulence.  It  is  a  measure  which 
depends  on  the  response  of  the  aircraft  to  the  turbulence,  a  response 
which  depends  upon  aircraft  weight,  altitude,  airspeed,  and  flight  charac¬ 
teristics.  The  derived  gust  velocity  attempts  to  compensate  for  the  air¬ 
craft  parameters  which  affect  its  response  by  providing  an  estimate  of  a 
standardized  wind  gust  that  would  have  produced  the  measured  aircraft 
response.  In  particular,  the  aircraft  response  (vertical  acceleration) 
is  modeled  as  proportional  to  airspeed  and  inversely  proportional  to 

aircraft  weight  (Lee,  1974) . 

1/3 

In  this  study,  e  was  chosen  as  the  principal  measurement  of 

turbulence.  The  peak  vertical  acceleration  values  were  used  to  provide 

supporting  data  about  the  expected  pilot  response.  Five  separate  pilot 

response  estimates  are  displayed  in  Figure  22.  The  intensity  descriptors 

recommended  by  MacCready  (1964)  on  the  basis  of  a  study  by  Zbrozek  (1961) 

are  given  as  well  as  the  descriptors  proposed  by  Zbrozek.  The  subjective 

observations  reported  by  Ehernberger  (1974)  are  also  presented  together 

with  the  conventional  turbulence  level  criteria  adopted  by  the  U.S.  Dept. 

of  Commerce  Subcommittee  on  Aviation  Meteorological  Services  in  1966. 

The  T-28  pilot  reports  for  the  22  July  1976  penetration  flights  are  also 

1/3 

included.  They  are  displayed  vs  e  because  the  eddy  dissipation  rate 
was  used  for  the  primary  measurement  of  turbulence  intensity.  Reports 
were  not  used  for  the  22  June  case  studies  because  the  pilot  responded 
only  with  "good"  or  "heavy"  to  characterize  the  turbulence. 

Zbrozek  used  rms  gust  loading  as  his  quantitative  measure  of  air¬ 
craft  response  to  compare  with  pilot  reaction.  Because  the  peak  values 
used  in  this  study  were  the  largest  in  each  30-second  sample  of  flight 
segment  used  to  characterize  turbulent  intensity  levels,  the  rms  vertical 
acceleration  value  was  estimated  as  one-half  the  reported  (2  standard 
deviations)  value.  Ehernberger  employed  the  peak-to-peak  acceleration 


change  measurement  to  quantify  aircraft  response.  His  peak-to-peak 

values  were  reduced  by  a  factor  of  1.3  for  plotting  on  Figure  22.  The 

reduction  factor  is  the  median  factor  relating  the  observed  peak-to-peak 

changes  in  the  undersampled  acceleration  data  to  the  observed  peak 

deviations  from  quiet  conditions  during  the  22  July  1976  penetration 

flights.  The  reduction  factor  is  only  an  approximation  since  rapid 

fluctuations  are  known  to  occur  between  samples. 

From  Figure  22,  it  is  evident  that  a  wide  range  of  word  descriptors 

1/3 

are  possible  for  a  given  e  value.  Sand  et  al .  (1974),  in  reporting 

T-28  thunderstorm  penetrations  for  the  1972  NHRE  field  season,  adopted 

the  bounds  proposed  by  MacCready.  They  used  severe  rather  than  heavy  to 
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characterize  the  region  5  <  e  <12  cm  s  .  We  have  adopted  the 

word  descriptors  used  by  the  T-28  pilot  during  the  22  July  penetration 

flights.  The  eddy  dissipation  rate  bounds  employed  were  obtained  by 

1 ' 3 

comparing  the  pilot  reports  to  the  simultaneously  recorded  e  '  values: 

5  <_  <  12  is  described  as  moderate,  12  £  <  15  is  described  as 

1/3  2/3  -1 

severe,  and  e  £  15  cm  s  is  described  as  extreme.  The  several 
1/3  2/3  -1 

occurrences  of  e  greater  than  15  cm  s  were  accompanied  by  pilot 
reports  that  the  aircraft  was  barely  controllable  (Heymsfield  et  al . , 
1978).  This  condition  is  best  described  as  extreme,  in  keeping  with  the 
accepted  definition  of  extreme,  rather  than  by  severe  or  even  by  moderate 
(conventional  criterion,  see  Figure  22) . 

The  severity  of  the  turbulence  experienced  by  an  aircraft  depends 
upon  the  airspeed,  weight,  and  flight  characteristics  of  the  aircraft. 

If  the  derived  gust  velocity  model  is  employed  for  the  analysis  of 
severity,  aircraft  with  the  same  lift  parameters  will  experience 
increased  severity  with  increased  airspeed  and  decreased  severity  with 
increased  weight  in  encounters  with  the  same  intensity  of  atmospheric 
turbulence  (derived  gust  velocity) .  Because  the  weight  of  many  of  the 
commercial  transports  is  an  order  of  magnitude  or  more  than  the  weight 
of  the  T-28  but  the  cruise  velocity  is  only  three  to  four  times  that  of 
the  T-28,  inflight  turbulence  severity  levels  should  be  at  least  one 
intensity  class  lower  than  reported  for  the  T-28.  Conversely,  lighter 
general  aviation  aircraft  will  find  the  turbulence  to  be  more  severe. 
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4.2  Penetration  Flight  Summary 

1/3  2/3  -1 

All  occurrences  of  e  greater  than  5  cm  s  were  identified  and 

located  on  radar  cell  and  cluster  summary  position  plots  for  each  of  the 

penetration  flights.  The  flights  were  subdivided  into  15-second  segments 

1/3 

(approximately  1.5  km  of  flight)  and  the  largest  e  value  for  each 
interval  was  used  to  characterize  the  turbulence  within  a  segment.  A 
tally  was  kept  of  all  the  turbulent  segments  (patches)  and  their  intensity. 
Figure  23  presents  cumulative  distributions  of  the  turbulent  patches 
observed  by  encounters  (1)  along  the  entire  flight  path,  (2)  within  the 
7  dBZ  echo  region  observed  by  the  S-band,  CP-2  radar  plus  a  4-km  buffer 
outside  the  echo  boundary,  and  (3)  within  the  20  dBZ  echo  region  plus 
4-km  buffer. 

The  plotted  locations  of  the  turbulent  patches  encountered  on  the 

second  aircraft  penetration  on  22  July  1976  are  depicted  in  Figures  11 

and  12.  All  the  turbulence  encounters  for  the  22  July  case  study  are 

plotted  in  Figures  13c  and  13d;  in  these  figures,  values  in  excess 

2/3  -1 

of  5  cm  s  (moderate)  are  coded  by  cross-hatching  and  values  in  excess 
2/3  -1 

of  12  cm  s  (severe)  are  indicated  by  solid  bars  across  the  aircraft 
track.  The  locations  of  one  (or  more)  encounters  with  severe  turbulence 
are  also  indicated  in  Figure  14. 

The  times  of  the  turbulent  path  encounters  were  scaled  from  the 
computer  generated  strip  chart  plots  provided  by  SDSMT.  Figure  21  is  a 
reproduction  of  several  of  the  variables  displayed  on  the  computer 
generated  plots.  Daily  summaries  were  also  provided  by  SDSMT.  Appendix 
B  contains  the  daily  summaries  for  each  case  study. 
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PERCENTAGE  OF  FLIGHT  TIME 


ve  distributions  of  turbulent  patches  of  1.5  km  extent  encountered  during  660  In 
penetration  flights 


5.  RADAR-AIRCRAFT  DATA  COMPARISON: 

A  TEST  OF  THE  HAZARD  DETECTION  HYPOTHESIS 


The  cell  approach  to  hazard  detection  presumed  the  existence  of  a 
turbulent  patch  within  a  small  radius  of  influence  around  a  significant 
cell  (or  cluster  of  cells).  The  underlying  physical  hypothesis  suggests 
that  the  radius  of  influence  should  approximate  half  a  cell  diameter. 

The  statistical  analysis  of  the  average  cell  size  for  Kansas  summer 
storms  reported  by  Crane  and  Hardy  (1981)  indicated  that  the  cells  have 
an  exponential  size  distribution  with  an  average  diameter  of  2.5  km. 

The  actual  radius  of  influence  to  be  applied  to  hazard  detection 
should  be  larger  than  half  a  cell  diameter  because  the  statistical 
uncertainty  in  the  radar  reflectivity  of  a  resolution  element  produces 
an  attendent  position  uncertainty  in  the  centroid  location  of  a  cell. 

The  reflectivity  cell  encloses  all  the  contiguous  radar  resolution 
elements  having  reflectivity  values  within  3  dB  of  the  largest  value  for 
a  cell.  If  a  cell  is  at  a  range  of  50  km,  it  will  contain  four  resolu¬ 
tion  elements  on  average  (1°  beamwidth  radar  system).  If  a  limited 
number  of  independent  samples  are  used  to  calculate  the  reflectivity 
value,  the  statistical  uncertainty  in  the  reflectivity  values  will  affect 
the  number  of  elements  included  in  the  cell  on  any  one  scan.  It  was 
estimated  that  the  cell  centroid  position  uncertainty  should  increase 
the  radius  of  influence  by  about  1  km  at  a  range  of  50  km. 

Additional  uncertainties  produced  by  cell  tracking  errors  and  by 
aircraft  tracking  errors  should  increase  the  radius  of  influence  by 
about  another  kilometer.  The  a  priori  estimate  of  the  expected  radius 
of  influence  of  a  cell  was  between  2  and  4  km.  Inter-radar  comparisons 
of  reflectivity  and  tangential  shear  locations  (see  Figure  12  for  an 
example)  show  tightly  clumped  significant  cells.  The  clumps  of  signifi¬ 
cant  cells  observed  by  the  different  radars  can  generally  be  circum¬ 
scribed  by  a  circle  with  less  than  a  1.5  km  radius.  The  statistical 
positioning  uncertainty  is  therefore  less  than  expected  a  priori, 
especially  for  comparisons  between  reflectivity  and  tangential  shear 
cells  detected  by  different  radars. 

The  statistical  measures  used  to  test  the  hazard  detection  hypothesis 
were  evaluated  as  functions  both  of  turbulence  level  and  of  radius  of 
influence.  The  cell  detection  hypothesis  presumed  the  existence  of 
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significant  turbulence  but  of  an  unspecified  level.  In  the  statistical 

1/3 

analysis,  detection  is  counted  if  the  e  value  was  above  the  threshold 
level,  but  a  false  alarm  was  tallied  if  the  value  was  below  the  threshold. 
The  false  alarm  rates  therefore  must  increase  as  the  threshold  level  is 
increased. 

5.1  Probability  of  Detection 

A  successful  detection  of  a  turbulent  patch  was  counted  if  one  or 
more  significant  cells  were  within  the  specified  radius  of  influence  of 
the  turbulent  patch.  The  probability  of  detection  (POD)  was  calculated 
by 

POD  =  (X/W)  •  100  (percent)  (1) 

where  X  is  the  number  of  detected  turbulent  patches  (successes)  and  W  is 

the  total  number  of  turbulent  patches  (trials). 

The  probabilities  of  detection  are  plotted  versus  turbulence 
1/3 

intensity  threshold  (e  )  for  reflectivity  cells  detected  using  the 

CP-2  radar  (Figure  24a),  for  tangential  shear  cells  detected  by  the 

C-band  radars  (Figure  24b),  and  for  tangential  shear  cells  detected  by 

the  X-band  radars  (Figure  24c).  The  results  are  for  a  6-km  radius  of 

influence.  The  variation  of  detection  probability  with  radius  of 
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influence  for  turbulence  intensities  in  excess  of  5  cm  s  is  displayed 

in  Figure  1.  Detailed  tabulations  of  the  detection  probabilities  by 

radar,  radius  of  influence,  and  threshold  level  are  included  in  Appendix  C. 

Three  curves  are  presented  in  Figure  24a  for  the  CP-2  radar,  one 

for  using  all  the  volume  cells,  the  second  for  volume  cells  with  an 

average  reflectivity  in  excess  of  40  dBZ  and  the  third  for  volume  cells 

with  an  average  reflectivity  in  excess  of  50  dBZ.  The  detection  proba- 

1/3 

bilities  (POD)  increase  with  increasing  threshold  level  for  5  <  e  < 

2/3  -1  ~ 

11  cm  s  They  decrease  with  increasing  average  reflectivity  values. 

Figure  1  displays  the  same  dependence  on  reflectivity  and  indicates  that 

the  larger  the  radius  of  influence,  the  higher  the  probability  of 

detection. 

The  data  presented  in  Figure  24  are  for  all  encounters  with  turbu¬ 
lence;  the  data  presented  in  Figure  1  are  for  encounters  within  regions 
containing  echoes  greater  than  20  dBZ  as  observed  with  the  CP-2  radar 


POD  {%) 


*00  <%J 


Figure  24c  Probability  of  detecting  a  turbulent  patch  within  6  kn  of  a 
significant  cell 
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plus  a  4-km  buffer  zone.  Appendix  C  lists  POD  values  for  all  encounters, 
encounters  within  the  7  dBZ  echo  region  plus  buffer  zone,  and  encounters 
within  the  20  dBZ  region  plus  buffer  zone.  The  size  of  the  buffer  zone 
was  based  on  the  radius  of  influence.  Assuming  that  a  cell  has  a  radius 
of  between  1  and  2  km,  the  closest  a  cell  can  be  to  the  edge  of  an  echo 
region  is  2  km.  The  buffer  zone  is  2  km  smaller  than  the  radius  of 
influence  and  extends  the  detection  region  to  the  maximum  possible  for 
cells  detected  within  a  region  with  echoes. 

The  statistics  for  turbulent  patch  detections  for  the  Doppler  radars 
are  presented  in  Figures  24b  and  24c.  The  results  for  the  CP-2  radar 
using  all  the  observed  cells  are  also  presented  in  these  figures  for 
comparison.  No  significant  differences  are  evident  between  the  results 
obtained  using  all  the  reflectivity  derived  significant  volume  cells  and 
the  results  obtained  using  the  single  Doppler  C-band  radars.  CP-3 
appeared  to  have  performed  better  than  CP-2  or  CP-4,  but  data  for  CP-3 
were  available  only  from  the  22  July  case  study.  The  POD  values  were 
significantly  lower  for  the  22  June  day  than  for  the  22  July  day  because 
the  aircraft  flight  included  two  cloud  penetrations  (Cl  and  C2,  see 

Appendix  B)  on  that  day  which  included  22  turbulent  patches.  Two  of  the 

1/3  2/3  -1 

within  cloud  patches  reached  an  e  value  in  excess  of  11  cm  s 

The  X-band  Doppler  radars  displayed  lower  POD  values  at  the  lower 
1/3 

e  threshold  levels  but  produced  results  equivalent  to  the  other  radars 
at  the  higher  levels.  The  lower  POD  values  are  attributed  to  attenuation. 
The  aircraft  penetration  flights  tended  to  be  along  the  upwind  flank  of 
the  storm  which  was  the  other  side  of  the  storm  from  the  X-band  radars 
(see  Figure  11).  The  intervening  rain  apparently  reduced  the  reflectivity 
values  sufficiently  to  prevent  detection  at  the  17  dBZ  reflectivity 
threshold  level  employed  in  the  data  processing. 

5.2  False  Alarm  Rate 

The  results  of  the  probability  of  detection  analysis  suggest  the 
use  of  a  large  radius  of  influence.  A  very  large  radius  would  defeat  the 
intent  of  the  cell  detection  approach  to  localizing  the  regions  of 
potential  hazard  and  should  create  a  high  false  alarm  rate.  False  alarms 
(failures)  were  defined  as  aircraft  penetrations  within  the  radius  of 
influence  of  a  significant  volume  cell  that  were  not  turbulent  at  or 
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above  the  expected  threshold  level  for  z  '  .  The  false  alarm  rate  (FAR) 
was  calculated  by: 

/Y  "  Y  \ 

FAR  =  I  -  1  •  100  (percent)  (2) 

where  Yq  is  the  total  number  of  significant  volume  cells  within  the 
specified  radius  of  influence  of  the  aircraft  track  and  Y  is  the  number 
of  volume  cells  within  the  specified  radius  of  influence  of  a  turbulent 
patch  (success). 

The  false  alarm  rates  for  each  of  the  radars  are  displayed  as  a 
function  of  threshold  level  in  Figure  25.  As  expected,  the  FAR  increases 
with  threshold  level.  However,  as  shown  in  Figure  2,  the  FAR  decreased 
with  increasing  radius  of  influence.  It  is  anticipated  that  the  FAR 
value  reaches  a  minimum  then  increases  as  the  radius  of  influence  con¬ 
tinues  to  increase. 

The  aircraft  penetration  flights  were  all  in  regions  of  new  and 
developing  convective  cells.  If  the  flights  had  been  through  the  older 
regions  containing  debris  cells  (see  Figure  4d;  see  also  Crane  and  Hardy, 
1981)  the  false  alarm  rate  could  have  been  higher. 

5.3  Critical  Success  Index 


The  optimum  hazard  detection  technique  would  maximize  the  detection 
probability  and  minimize  the  false  alarm  rate.  Critical  success  indicies 
are  often  used  to  combine  both  detection  (success)  and  false  alarm 
(failure)  data  in  a  single  figure  of  merit  for  forecast  or  detection 
schemes.  A  critical  success  index  was  adopted  for  this  study  which 
focused  on  the  behavior  of  the  cell  approach  to  hazard  detections.  It 
counted  cells  that  were  associated  (within  a  radius  of  influence)  with  a 
turbulent  patch  as  a  success  and  a  turbulent  patch  that  was  not  detected 
as  a  failure.  The  resultant  critical  success  index  (CSI)  is  given  by: 


CSI 


Y  +  W  -  X 


100  (percent) 


(3) 


v°  / 

where  the  parameters  are  defined  for  equations  (1)  and  (2). 

The  CSI  values  are  plotted  as  a  function  of  turbulence  threshold  in 
Figure  26  and  as  a  function  of  radius  of  influence  in  Figure  3.  The 
results  indicate  a  decrease  in  the  performance  of  the  hazard  detection 
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Figure  25b  Probability  of  a  false  alarm  (no  turbulence)  along  the  entire 
segment  of  aircraft,  track  within  6  km  of  a  significant  cell 
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Figure  25c  Probability  of  a  false  alarm  (no  turbulence")  along  the  entire 
segment  of  aircraft  track  within  6  km  of  a  significant  cell 
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Figure  26c  Critical  success  index 
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hypotheses  (CSI  value)  as  the  turbulence  threshold  is  increased  but  an 
increase  in  CSI  as  the  radius  of  influence  is  increased. 

The  results  for  reflectivity  data  alone  show  that  the  CSI  is  a 
1/3  2/3  -1 

maximum  at  the  lowest  e  value  (5  cm  s  )  when  all  reflectivity 

1/3  2/3  -1 

levels  are  used  but  is  a  maximum  at  e  'v  10  cm  s  when  only  the 
higher  reflectivity  values  are  employed.  The  use  of  the  higher  reflec¬ 
tivity  values  alone  does  not  improve  the  performance  of  the  algorithm 
relative  to  the  results  obtained  using  all  the  data. 

The  use  of  Doppler  data  appears  to  improve  the  hazard  detection 
results  (CSI  values)  in  the  8  <_  <_  12  cm^^  s  *  range  relative  to  the 

results  obtained  using  reflectivity  data  alone.  At  lower  thresholds, 
the  C-band  Doppler  data  and  the  reflectivity  data  alone  produce  comparable 
results  (Figure  3) . 

The  results  obtained  for  a  6-km  radius  of  influence  and  a  >_ 
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5  cm  s  threshold  are  excellent  as  illustrated  in  Table  1.  The 
curves  in  Figures  1  through  3  suggest  even  better  results  for  larger 
radii  of  influence.  The  mid-level  flow  field  depicted  in  Figure  15b 
indicates  that  a  distance  as  large  as  10  km  may  be  required  to  embrace 
regions  of  blocked  flow  upwind  of  a  strong  updraft  especially  as  the 
region  develops  and  few  tracers  are  available  to  provide  radar  detection. 

The  radius  of  influence  analysis  was  not  extended  to  distances 
beyond  6  km  for  several  reasons.  The  intent  of  the  analysis  was  to  show 
that  the  cell  approach  to  hazard  detection  was  viable  and  a  radius  of 
influence  did  exist.  The  intent  of  the  analysis  was  satisfied  by  the 
results.  The  selection  of  an  optimum  radius  of  influence  depends  upon 
the  environment  through  which  penetrations  are  to  be  attempted.  The 
values  may  be  different  for  flights  through  the  debris  downwind  of  a 
cell  from  flights  through  the  regions  of  new  growth.  The  former  regions 
should  be  less  turbulent  and  are  therefore  of  more  interest  to  air 
traffic  control.  Measurements  should  be  made  under  the  conditions  for 
which  the  system  should  operate  before  an  optimum  set  of  parameters  can 
be  established.  More  than  two  sets  of  observations  are  required  to 
establish  the  optimum  set  of  parameters. 
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The  aircraft  penetrations  flown  on  22  July  were  used  to  compare  the 

observed  turbulence  levels  with  the  significant  volume  cell  parameters. 

1/3 

The  highest  level  (e  )  patch  within  the  specified  radii  of  influence 

are  listed  for  each  cell  observed  within  that  radius  of  the  aircraft 

track  in  Appendix  A.  Figure  27  presents  the  observed  relationships 

1/3 

between  volume  cell  average  reflectivity  and  e  ;  Figure  28  presents 

the  observed  relationship  between  the  averaged  magnitude  of  the  tangential 
1/3 

shear  and  e  .  No  correlation  is  evident  in  either  figure.  If  anything, 
the  highest  observed  value  was  associated  with  the  lowest  observed 

tangential  shear  values.  • 

A  lack  of  correlation  between  derived  gust  velocities  and  colocated 
reflectivity,  reflectivity  gradient,  Doppler  spread,  or  Doppler  velocity 
shear  values  was  reported  by  Lee  (1977).  He  associated  the  radar  obser¬ 
vations  and  turbulence  encounters  over  distances  less  than  2  km.  A 
similar  lack  of  correlation  is  reported  here  where  the  associations  are 
to  local  maxima  in  the  reflectivity  or  shear  fields  at  distances  up  to 
6  km. 

The  results  obtained  to  date  indicate  that  cell  age  is  more  important 
than  reflectivity  or  tangential  shear  in  establishing  a  correlation 
between  e^"5  and  the  cell  parameters.  The  several  cases  of  severe 
turbulence  depicted  in  Figure  13c  and  13d  all  occur  outside  regions  of 
high  reflectivity  but  are  associated  with  regions  that  later  produced 
significant  cells  and  clusters. 

Significant  regions  of  turbulence  are  most  active  just  before  or  at 
the  time  of  first  detection.  They  cannot  be  detected,  processed,  and 
used  to  forecast  hazard  location  10  to  20  minutes  in  the  future  because 
they  did  not  exist  for  a  long  enough  time  to  be  tracked  before  the 
region  becomes  hazardous.  No  radar  based  system  of  observations  will 
solve  this  basic  hazard  detection  problem.  Only  forecasts  of  probable 
locations  of  new  cell  (or  cluster)  development  can  provide  a  solution  to 
the  hazard  detection  problem.  The  regular  organization  of  new  cluster 
development  (Figure  14)  indicates  that  the  larger-than-cloud  scale 
organization  controls  the  development  of  new  cells  and  that  an  early 
recognition  of  the  pattern  of  new  cell  development  could  be  useful  in 
providing  forecasts. 
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Figure  28  Observed  relationship  between  eddy  dissipation  rate  and 
volume  cell  tangential  shear,  22  July  1976 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  documents  the  development  and  testing  of  a  viable  air¬ 
craft  hazard  detection  scheme  which  operates  as  well  employing  conventional 
radar  data  as  it  does  employing  Doppler  radar  data.  The  detection  scheme 
was  evaluated  using  aircraft  penetrations  into  the  most  turbulent  regions 
of  hail  storms.  These  regions  are  of  interest  when  verifying  that 
turbulence  can  be  detected  but  are  not  optimum  for  evaluating  the  detec¬ 
tion  technique  under  realistic  operating  conditions.  If,  in  the  latter 
case,  higher  false  alarm  rates  can  be  expected,  more  development  will  be 
necessary  to  utilize  more  effectively  the  available  information  before 
the  hazard  detection  procedure  can  be  incorporated  into  an  air  traffic 
control  environment. 

Continued  algorithm  refinement  is  necessary  and  should  be  pursued. 

The  horizontal  flow  field  maps  (Figure  15)  showed  that  single  Doppler 
radars  can  be  successfully  employed  to  locate  the  cores  of  the  strongest 
(non-entraining)  updrafts.  Cluster  age  was  also  an  important  indicator 
of  severity.  It  is  an  indicator  which  will  be  important  in  separating 
cells  or  clusters  formed  in  the  debris  region  downwind  of  an  active 
cluster  from  the  active  and  potentially  hazardous  growth  regions.  The 
progressive  cluster  development  pattern  can  be  used  to  separate  automati¬ 
cally  new  growth  and  debris  regions. 

Continued  flight  testing  should  proceed  hand-in-hand  with  the  con¬ 
tinued  algorithm  development.  Of  particular  importance  is  the  develop¬ 
ment  of  a  real-time  significant  cell,  cluster,  and  potential  hazard 
display  capability.  The  real-time  display  provides  the  only  means  to 
evaluate  successfully  the  hazard  detection  scheme.  Aircraft  must  be 
vectored  into  apparently  hazard-free  areas  to  evaluate  the  probability 
of  detection  and  into  apparently  hazardous  areas  to  evaluate  the  proba¬ 
bility  of  a  false  alarm.  Only  with  positive  flight  control  will  enough 
data  be  amassed  to  optimize  threshold  criteria,  the  radius  of  influence, 
and  algorithm  improvements  yet  to  be  made. 

A  major  conclusion  of  the  work  reported  here  is  that  the  severity 
of  the  turbulence  was  inversely  related  to  the  age  of  the  cluster;  the 
younger  the  cluster,  the  more  turbulent.  Unfortunately,  a  detection 
scheme  based  only  on  observations  requires  several  observations  to 
establish  the  existence  and  motion  of  a  cluster.  The  danger  is  highest 
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before  the  hazard  can  be  detected.  A  short  range  forecast  of  the  probable 
locations  of  new  cell  development  is  required  to  complete  the  development 
of  the  hazard  detection  scheme. 

The  recommended  program  to  develop  the  automatic  hazard  detection 
system  for  use  by  the  FAA  is: 

1)  refine  the  cell  and  cluster  detection  algorithms  to 
reduce  the  effect  of  the  statistical  uncertainties  in 
the  estimation  of  velocity  perturbations; 

2)  refine  the  cell  significance  algorithms  using  cell  age, 
location  in  the  development  pattern,  apparent  vertical 
transport  of  horizontal  momentum  and  similar  parameters 
to  provide  an  indication  of  the  intensity  of  the  tur¬ 
bulence; 

3)  develop  a  graphical  display  of  the  cell,  cluster,  and 
contour  data  that  can  be  provided  both  to  meteorologists 
for  the  development  of  short  range  forecasts  and  to 
controllers  for  use  in  vectoring  test  aircraft  into 
regions  expected  to  be  hazardous  or  free  of  hazard  to 
evaluate  the  algorithms; 

4)  develop  a  real  time  version  of  the  program  for  use  in 
the  continued  evaluation  of  the  algorithms; 

5)  using  the  real  time  system  and  the  controller  displays, 
conduct  flight  tests  to  evaluate  the  refined  velocity 
perturbation  detection  and  turbulence  intensity  esti¬ 
mation  algorithms; 

6)  analyze  the  new  cell  and  cluster  development  patterns 
to  isolate  the  propagating  disturbances  that  trigger 
the  convection; 

7)  automate  the  display  of  the  propagating  disturbances 
to  provide  forecasters  and  controllers  with  the  tools 
necessary  for  the  short  range  forecast  of  hazardous 
areas; 

8)  evaluate  the  short  range  hazardous  region  forecast 
procedures  in  several  different  climate  regions;  and 
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9)  based  on  the  short  range  hazard  forecast  evaluations, 
automate  the  successful  features  to  provide  timely 
forecasts  to  the  controllers. 

The  recommended  program  includes  tasks  that  can  be  readily  defined 
at  the  present  time  (1  through  4  and  6),  tasks  to  be  completed  when 
refined  algorithms  are  available  (5  and  8),  and  tasks  whose  definition 
depends  on  the  results  of  future  work  (7  and  9). 
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19  January  1978 


MEMORANDUM  FOR  THE  RECORD 


FROM: 


P.  L.  Smith,  Jr. 


SUBJECT:  Notes  Based  on  Perusal  of  T-28  Data  Plots  for  22  June  1976 
(Flight  167) 


All  numerical  values  given  are  approximate  values  read  from  the 
computer-generated  plots.  Detailed  listings  should  be  consulted  if  more 
precise  values  are  desired. 


la  S  lb  Penetration  1  (l6:02:0U  to  16:09:27) 

Altitude  range:  6.1  to  6.5  tan 
Temperature  range:  -IOC  to  -16.5C 
Updraft  characteristics: 

Maximum  updraft  speed  about  lU  m  sec-1  near  16:05:36. 

Width*  of  main  updraft  about  5  tan. 

Gradient  of  0e  across  main  updraft  8K  (partly  related  to 
altitude  change? ) . 

Maximum  downdraft  speed  about  10  m  sec-1  (two  locations). 

Width*  of  main  downdrafts  3  -  U  km. 

Hydrometeor  characteristics: 

Cloud  liquid  water  concentration  up  to  about  0.9  g  m~3  except 
for  suspicious  looking  peaks  between  1603  and  I60U.  Little 
cloud  liquid  water  in  the  main  updraft  (perhaps  because  of 
depletion  by  the  hail?).  Cloud  LWC  values  up  to  about  0.7  g  m” 
appear  in  some  downdraft  regions. 

Hail  present  in  substantial  concentrations  in  the  main  updraft 
at  1605-1606  as  well  in  the  downdraft  region  around  1607.  Hail 
recorded  continuously  for  more  than  40  km  (l6oU-l6ll+) .  Hail 


*  Extent  of  the  draft  region  along  the  flight  path,  which  is  probably 
smaller  than  its  true  width. 
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mass  concentrations  of  more  than  10  g  m~3  are  suspicious*,  as 
are  reflectivity  factors  up  to  73-7**  dBz.  Hail  concentrations 
exceeding  1  g  m-3  show  up  outside  the  cloud  (16:09:27  -  16-.11+); 
pilot  comments  should  be  checked  to  see  whether  there  is  any 
corresponding  annotation. 

Turbulence:  Maximum  values  about  10  cm2/3  sec-1,  occuring  at 
edges  of  main  updraft . 

Acceleration:  Less  than  0.5  g  either  direction.  Maximum  values 
occur  in  downdrafts . 

Cl  §  C2  In-Cloud  Period  Between  16:13:30  and  16:25:00  or  Later 

This  period  which  had  a  peak  updraft  exceeding  20  m  sec-1  and  cloud 
liquid  water  concentrations  up  to  1.5  g  m-3  was  not  logged  as  a 
penetration.  During  this  period,  there  are  several  places  where 
the  reverse  flow  temperature  (RFT)  probe  trace  and  the  Johnson- 
Williams  (JW)  liquid  water  concentration  trace  look  suspisciously 
similar.  See  for  example  l6 :1** :25-l6:15 :**0  or  16: 18:00-16:18:30  or 
16:19:25-16:20:00  or  l6: 21: 55-16:23: 00  or  16:2** : 10-16:25 : 05.  There 
is  little  suggestion  of  similar  behavior  during  Penetration  1.  An 
acceleration  peak  of  +0.7  g  occurred  at  a  "hole  in  an  updraft"  just 
after  16:19:1*5. 

3a  5  3b  Penetration  2  (l6: 25: 1*1-16:39: 29) 

Altitude  range:  5.2  to  6.7  km 
Temperature  range:  — Uc  to  -17C 
Updraft  characteristics: 

Maximum  updraft  speed  about  16  m  sec-1  (several  places;  probably 
underestimated  due  to  effects  of  aircraft  icing). 

Width  of  main  updraft  about  8  km;  actually  about  30  km  of 
essentially  sustained  upward  motion  from  16:27:1*5  to  16:32:1*5. 
During  the  5  min  crossing  of  the  large  (30  km)  updraft  area, 
the  aircraft  gained  about  1.35  km  in  altitude.  The  peak 
updraft  speed  was  only  about  15  m  sec-1 ,  but  the  average  for 


*  But  note  that  the  hail  concentrations  plotted  are  for  1-sec  samples; 

they  should  be  averaged  over  several  seconds  to  get  representative  values. 
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the  entire  30  km  was  4.5  m  sec-1.  Cloud  liquid  water 
concentrations  were  so  large  and  extensive  that  the  aircraft 
picked  up  enough  ice  to  require  descending  below  the  melting 
level  to  thaw  out. 

Maximum  downdraft  speed  about  22  m  sec-1  near  16:34:05 
(probably  overestimated). 

Width  of  main  downdraft  areas  about  3  km  (indications  of 
broader  downdrafts  after  1637  believed  due  to  degradation 
of  aircraft  performance  because  of  structural  icing). 

Hydrometeor  characteristics: 

Cloud  liquid  water  concentrations  reached  about  2  g  m-3;  values 
of  0.8  g  m-3  or  more  were  indicated  across  the  full  30  km  width 
of  the  major  updraft  area,  and  indications  of  cloud  LWC  above 
zero  were  maintained  for  more  than  60  km.  Cloud  liquid  water 
concentrations  of  more  than  1  g  m-3  appear  in  downdrafts  in 
several  places. 

Hail  concentrations  up  to  about  9.6  g  m-3,  mostly  in  updraft 
regions.  Width  of  main  hail  zone  about  13  km  (2:10);  inter¬ 
mittent  hail  encounters  over  more  than  40  km  (l6:26:45-l6:33:45) . 
Significant  amounts  of  hail  coexist  with  more  than  1  g  m-3  of 
cloud  liquid  in  the  major  updraft  area.  Reflectivity  factors 
up  to  68  dBz  are  more  plausible. 

Turbulence:  Major  jolt  of  28  cm2/3  sec-1  at  16:27:40;  two  other 
peaks  in  excess  of  15.  The  major  jolt  occurred  at  a  sharp  -14  m  sec-1 
to  +10  m  sec-1  downdraft /updraft  boundary;  the  secondary  peaks  were 
in  updraft  regions. 

Acceleration:  Peaks  in  both  directions  of  nearly  1  g  indicated 
around  1628;  several  peaks  of  more  than  +_  0.5  g. 

The  correspondence  between  the  RFT  and  JW  traces  is  again  remarkable; 
see  for  example,  16:31:00-16:35:00  or  later.  The  temperature  excursions 
at  16:33-16:34  are  especially  implausible  and  the  strong  correlation  between 
the  RFT  and  J V/  traces  during  that  period  is  highly  suspect. 

The  period  after  about  16:37  was  the  descent  to  allow  for  deicing 
the  aircraft.  During  this  period,  the  plot  shows  downdrafts  which  on 
occasion  exceed  15  m  sec-1 .  Our  updraft  formula  is  supposed  to  correct 
for  aircraft-induced  effects,  and  thus  there  should  be  essentially  no  up 
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or  downdraft  in  this  region  of  the  plot.  This  suggests  that  there  may 
be  some  difficulty  with  the  calculation.  One  factor  is  doubtless  that 
the  amount  of  ice  accumulated  degraded  the  aircraft  performance  so  that 
with  the  indicated  airspeed  and  manifold  pressure  at  normal  level  flight 
conditions,  the  T-28  would  actually  be  sinking;  our  computations  would 
treat  the  sink  rate  as  a  downdraft.  In  support  of  this  argument,  the 
departures  from  zero  indicated  vertical  air  motion  diminished  markedly 
within  a  couple  of  minutes  after  the  aircraft  descended  below  the  OC  level. 
It  is  also  worth  noting  that  the  K-probe  was  carried  on  this  flight,  adding 
to  the  effects  of  icing  on  aircraft  performance. 

4a  4  4b  Penetration  3  (l6: 48:20-16: 36: 51 ) 


Altitude  range:  5.0  to  5-6  km 
Temperature  range;  -3C  to  -8C 
Updraft  characteristics: 

Maximum  updraft  speed  about  16  m  sec"1  at  16:50:20. 

Width  of  main  updraft  about  4  km. 

Maximum  downdraft  speed  about  8  m  sec-1. 

Hydrometeor  charcteristics : 

Cloud  liquid  water  concentrations  up  to  l.T  g  m-3;  peak  value 
observed  on  the  edge  of  a  downdraft  region. 

No  hail  recorded  on  this  penetration.  Pilot  comments  that 
aircraft  was  in  and  out  of  cloud,  and  that  no  hail  was 
encountered. 

Turbulence:  Exceeded  15  cm2/3  sec-1  at  16: 50:40  in  a  downdraft 
region.  Another  indicated  peak  at  about  16:49:10  is  suspicious 
because  its  appearance  and  magnitude  (near  4o)  suggest  an  error 
of  some  kind. 


Acceleration:  Maximum  +0.6  g  near  16:50:45  at  downdraft /updraft 
boundary . 
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The  temperature  probe  still  shows  evidence  of  following  the  JW  LWC 
trace  (see  l6:U9:*+5  and  later). 
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MEMORANDUM  FOR  THE  RECORD 
FROM:  P.  L.  Smith,  Jr. 


SUBJECT:  Notes  Based  on  Perusal  of  T-28  Data  Plot  and  Related  Data 
for  22  July  1976  (Flight  183) 


Most  numerical  values  given  are  approximate  values  read  from  the 
computer-generated  plots.  Detailed  listings  should  be  consulted  if 
more  precise  values  are  desired. 

1  Penetration  1  (16:30:31  to  l6 : 3U : ^3 ) 


Altitude  range:  6.1  to  6.7  km 
Temperature  range:  -10. 5C  to  -15. 5C 
Updraft  characteristics: 

Maximum  updraft  speed  21  m  sec-1  near  16:32:55. 

Width*  of  updraft  about  6  km. 

Range  of  8£  across  updraft  about  10  K. 

Maximum  downdraft  speed  6  ra  sec-1  at  16:33:1**. 

Width*  of  main  downdraft  about  4  km. 

Hydrometeor  characteristics: 

Cloud  liquid  vater  concentrations  up  to  2.7  g  m~3  in  updraft. 

Foil  impactor  data  show  significant  numbers  of  particles  larger 
.  than  3  mm  from  about  16:32:23  to  16:33:31  (approximately  7  km, 
mostly  in  updraft).  Counts  start  about  15  seconds  later  and 
continue  about  **5  seconds  longer  than  hail  spectrometer  counts, 
suggesting  the  possibility  of  a  time  shift  in  the  foil  analysis 
(although  differences  in  sampling  volumes  and  size  ranges  are 
other  possible  explanations.) 


•Extent  of  the  draft  region  along  the  flight  path,  which  is  probably  smaller 
■  than  its  true  width. 
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Maximum  recorded  particle  size  nearly  1.5  cm 
around  l6:32:U5  (in  updraft)  and  again  near 
16:33:21  {in  downdraft). 

Maximum  number  concentration  ( D  >  3  mm)  10-20/m3 
in  updraft  from  about  16:32:30  to  16:33:10. 

Maximum  indicated  mass  concentration*  more  than 
3  g/m3  from  roughly  l6:32:b0  to  16:32:50. 

Hail  spectrometer  active  from  16:32:10  to  l6:32:U3 
(approximately  3  km,  mostly  on  the  weak  side  of  the 
updraft ) . 

Maximum  hail  mass  concentration  0.7  g/m3  at 
16:32:33  in  updraft  with  about  0.3  g/m3  cloud 
LWC . 

Maximum  observed  hailstone  size  about  1.0  cm 
(Category  No.  6;  two  places). 

Maximum  number  concentration  7/m3,  also  at  16:32:33. 

Possible  indications  of  hail  depletion  of  cloud 
water  in  updraft  around  l6: 32: 10-16:32: 35- 

Reflectivity  factors  to  about  50  dBz. 

Turbulence:  Maximum  intensity  11  cm2/3  sec-1  in  the  updraft. 

Acceleration:  Peaks  about  +0.5  g  in  and  on  edges  of  the  updraft; 
a  recorded  -0.7  g  was  not  plotted. 

There  is  some  indication  that  the  reverse  flow  temperature  (RFT)  trace 
follows  the  Johnson-Williams  liquid  water  concentration  (LWC)  trace  from 
about  l6 : 32 : 20  onwards . 

2_  Penetration  2  (16:33:56  to  l6:U5:^8) 

Altitude  range:  6.2  to  6.9  km 


Foil-derived  mass  concentrations  quoted  are  12-sec  averages  for  particles 
larger  than  3  mm;  contribution  of  1-3  mm  particles  on  this  flight  is  always 
less  than  0.26  g/m3. 
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Temperature  range:  -12. 5C  to  -18.5C 
Updraft  characteristics: 

Maximum  updraft  speed  of  33  m  sec-1  at  l6:U0:U0. 

Width  of  main  updraft  about  1*  km. 

Range  of  0e  across  main  updraft  about  5  K. 

Maximum  downdraft  speed  17  m  sec-1  at  l6:Ul:30. 

Width  of  main  downdraft  about  5  km. 

Hydrometeor  characteristics: 

Cloud  liquid  water  concentrations  up  to  2.3  g  m-3 
predominantly  in  main  updraft;  LWC  values  to  O.U  g  m“3 
in  some  downdraft  regions. 

Foil  data  show  significant  numbers  of  particles  larger 
than  3  mm  from  about  l6:U0:12  to  l6: 1*1: 1*1*  (approximately 
9  km  across  main  updraft  and  downdraft). 

Maximum  recorded  particle  size  about  0.8  cm  near 
peak  of  main  updraft  around  16:1*0:38. 

Maximum  number  concentration  ( D  >  3  mm)  about 
l6/m3  in  main  downdraft  near  16:1*1:20. 

Maximum  indicated  mass  concentration  about  0.5  g/m3 
in  main  updraft  near  16:1*0:30  and  again  in  main  down- 
draft  around  l6:Ul:20. 

Hail  spectrometer  active  between  16:1*0:28  and  l6: 1*1:39 
(about  7  km),  with  a  few  sporadic  counts  elsewhere. 

Width  of  main  hail  zone  about  5  km,  on  edge  of  main 
updraft  and  primarily  in  main  downdraft. 

Maximum  hail  mass  concentration  1.9  g/m3  at 
16:1*1:20  in  main  downdraft. 

Maximum  observed  hailstone  size  nominally  1.25  cm 
(Category  8;  two  places). 
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Maximum  number  concentration  32/m3  at  16:1*1:21. 

Suggestion  of  hail  depletion  of  cloud  liquid  water 
in  region  of  main  updraft  after  16:1*0:1*5. 

Reflectivity  factors  to  55  dBz  are  plausible. 

Turbulence:  Maximum  intensity  11  cn2^3  sec-1  in  main  updraft. 

Acceleration:  Peak  +0.8  g  at  16:1*0:35  in  main  updraft. 

There  is  only  slight  indication  that  the  RFT  trace  follows  the  LWC 
trace  on  this  pass. 

Penetration  3  (l6;l*9:25  to  16:53:1*3) 

Altitude  range:  6.1  to  6.7  km 

Temperature  range:  -11C  to  -17C 

Updraft  characteristics: 

Maximum  updraft  speed  12  ra  sec-1  at  16:51:22  (not  in  main 
updraft ) . 

Width  of  main  updraft  about  5  km. 

Maximum  downdraft  speed  7  m  sec”1  around  16:51;  10  m  sec”1 
downdraft  20  seconds  (2  km)  before  cloud  entry. 

Downdraft  widths  of  the  order  of  2-3  km;  no  major  downdraft 
encountered  during  in-cloud  period. 

Hydrometeor  characteristics: 

Cloud  liquid  water  concentrations  to  2.1*  g  m"3  mostly  in 
updrafts,  but  values  to  near  1  g  m"3  occur  in  downdraft 
areas  as  well. 

Foil  irapactor  data  not  reduced;  few  particles  indicated. 

No  hail  encountered  on  this  penetration. 
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Turbulence:  Maximum  intensity  l6  cm2/3  sec-1  at  16:51:07  on 
edge  of  downdraft. 

Acceleration:  Peak  about  -0.5  g  at  16:52:35  in  main 

updraft  region;  recorded  peaks  of  +0.8,  -0.7  g  were 
not  plotted. 

From  16:50:30  the  RFT  trace  tends  to  follow  the  JW  LWC  trace  again, 
making  the  temperature  data  suspect. 

Penetration  U  ( 16:55:3**  to  17:01:59) 

Altitude  range:  6.1*  to  7.1  km 

Temperature  range:  -13. 5C  to  -21C 

Updraft  characteristics: 

Maximum  updraft  speed  25  m  sec-1  at  l6: 57:10. 

Width  of  main  updraft  about  5  km. 

Range  of  0e  across  main  updraft  about  7  K. 

Maximum  downdraft  speed  13  ra  sec-1  at  16:58:12. 

Width  of  main  downdraft  about  U  (7?)  km. 

Hydrometeor  characteristics: 

Cloud  liquid  water  concentrations  to  2  g  m-3,  essentially 
all  in  main  updraft . 

Foil  data  show  significant  numbers  of  particles  larger  than 
3  ram  from  about  16:57:30  to  16:59:02  (approximatley  9  km, 
predominantly  in  main  downdraft  region);  counts  start  about 
12  sec  later  than  hail  spectrometer  counts. 

Maximum  particle  size  recorded  about  1.1  cm  in 
downdraft  near  16:57:1*0  and  l6: 58:01*. 


Maximum  number  concentration  (D  >  3  mm)  20-23/m3 
in  main  downdraft  area  from  about  16:57:50  to  16:58:25. 
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Maximum  indicated  mass  concentration  about  1.8  g/ra3 
in  main  downdraft  near  16:58:00. 

Hail  spectrometer  active  between  l6: 57:18  and  16:58:5^ 
(approximately  9  km,  almost  entirely  in  main  downdraft 
region),  with  a  few  counts  elsewhere.  Width  of  main 
hail  zone  about  7.5  km. 

Maximum  hail  mass  concentration  more  than  10  g/m3 
at  16:57 :Ul  in  downdraft ;  value  of  W  essentially 
remains  above  1  g/m3  from  16:57:28  to  16:58:11. 

Maximum  observed  hailstone  size  nominally  U.U  cm 
(Category  No.  lU),  also  at  l6: 57:^1. 

Maximum  number  concentration  25/m3  at  16:57:55  in 
downdraft . 

Large  hail  mass  concentrations  and  reflectivity  factors 
up  to  72  dBz  are  suspicious.** 

Turbulence:  Maximum  intensity  12  cm2/3  sec-1  at  edge  of  main 
updraft . 

Acceleration:  Peaks  about  +0.5  g  in  and  at  edge  of  main  updraft; 
peak  of  +0.9  g  was  not  plotted. 

Little  indication  of  correlation  between  RFT  and  JW  LWC  traces. 

5_  Penetration  5  (17:0*»:U3  to  17:10:02) 

Altitude  range:  5.8  to  6.1  km 
Temperature  range:  -9C  to  -15- 5c 


**But  note  that  the  concentrations  plotted  are  for  1-sec  samples;  they 
should  be  averaged  over  several  seconds  to  get  representative  values. 


B-12 


Memo  For  The  Record  (PLS) 
7  February  1978 
Page  -  7  - 


Updraft  characteristics: 

Maximum  updraft  speed  8  m  sec-1  near  17:09. 

Updraft  areas  no  more  than  2  km  wide;  no  major  updraft 
encountered. 

Maximum  downdraft  speed  11  m  sec-1  near  17:07. 

Widths  of  main  downdraft  areas  3-**  km. 

Hydrometeor  characteristics: 

Very  little  cloud  liquid  water.  Maximum  cloud  LWC 
0.6  g  m“3;  some  smaller  values  in  downdraft  areas. 

Foil  data  show  significant  numbers  of  particles  larger 
than  3  mm  from  about  17:06:U3  to  17:09:03  (approximately 
lU  km  in  predominantly  downdraft  area) ;  counts  start 
about  10  sec  later  and  continue  about  10  sec  longer 
than  hail  spectrometer  counts,  suggesting  the  possibility 
of  a  time  shift  in  the  foil  analysis. 

Maximum  particle  size  recorded  about  1.0  cm  near 
17:08:13. 

Maximum  number  concentration  (D  >  3  mm)  10-lU/m3 
from  about  17:08:07  to  17:08:51. 

Maximum  indicated  mass  concentration  about  1.3  g/m3 
from  17:08:35  to  17:08:U7  (with  a  secondary  max 
0. 9-1.0  g/m3  around  17:08:03-17:08:15). 

All  these  maxima  occurred  in  downdrafts. 

Hail  spectrometer  active  from  17:06:31*  to  17:08:52 
(approximately  lU  km,  essentially  all  in  downdrafts). 

Maximum  hail  mass  concentration  about  5  g/m3  at 
17:08:lU;  value  of  W  remained  above  1  g/m3  from 
17:08:10  to  17:08:3>». 
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Maximum  observed  hailstone  size  nominally  3.**  cm 
(Category  No.  13)  at  17:08:03  and  17:08:10. 

Maximum  number  concentration  27 /m3  at  17:08:15. 

Reflectivity  factors  to  67  dBz. 

Turbulence:  Maximum  value  8  cm2/3  sec”1  in  downdraft. 

Acceleration:  PeaX  +0.5  g  in  downdraft. 

(Note  to  Bob  Crane:  Slight  baseline  shift  occured  on  the  plotter 
at  about  17 : 08 : 10 . ) 

6  Penetration  6  (17:13:41  to  17:22:12) 


Altitude  range:  5.6  to  6.4  km 

Temperature  range:  -9C  to  -15C 

Updraft  characteristics: 

Maximum  updraft  speed  21  m  sec”1  at  17 :  1*» :  ^*5 . 

Width  of  main  updraft  about  5  km. 

Maximum  downdraft  speed  20  ra  sec-1  at  17:15:30. 

Width  of  main  downdraft  2  (or  4?)  km. 

Hydrometeor  characteristics: 

Cloud  liquid  water  concentrations  up  to  2.3  g  m”3  in 
updrafts ;  values  up  to  more  than  1  g  m"3  in  downdrafts . 
Cloud  LWC  above  zero  for  more  than  2  min  (12  km). 

Foil  impactor  Jammed  early  in  penetration,  apparently 
because  of  ice  accumulated  during  previous  penetrations; 
data  not  reduced. 

Hail  spectrometer  intermittently  active  from  17:14:27 
to  17:22:04  (about  45  km),  mainly  in  downdraft  areas. 
Substantial  amounts  of  hail  found  only  in  major  down- 
draft  areas  between  17:15:03  and  17:16:24;  main  hail 
zone3  about  2  km  wide.  Sporadic  indications  of  small 
hail  elsewhere. 
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Maximum  hail  mass  concentration  2.5  g/m 3  at 
17:16:02  on  edge  of  downdraft. 

Maximum  observed  hailstone  size  about  2.0  cm 
(Category  No.  11 )  also  at  17:16:02. 

Maximum  number  concentration  27/m3  at  same 
location. 

Hail  and  nearly  1  g  m“3  of  cloud  liquid  water 
appear  together  in  a  downdraft  region  after 
17 : l6 . 

Reflectivity  factors  to  55  dBz  are  plausible. 

Turbulence:  Maximum  intensity  l6  cm2^3  sec-1  at  peak  downdraft 
at  17:15:30. 

Acceleration:  Peaks  +0.9,  -0.7  g  in  downdraft  around  17:15:30; 
a  -0.8  g  peak  was  recorded  but  not  plotted. 

The  RFT  trace  conspicuously  follows  the  JW  LWC  trace  after  17:1^:10. 
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APPENDIX  C 


STATISTICAL  HAZARD  DETECTION  SUMMARIES 
BY  RADAR  AND  RADIUS  OF  INFLUENCE 


Radar: 


7  d BZ  20  dBZ 


C-10 


Radar:  CP-3,  Doppler  Radius:  2  km 


Radar:  CP-3,  Doppler  Radius:  4  km 


Radar:  CP-3,  Doppler  Radius:  6  km 


Radar:  CP-4,  Doppler  Radius:  2  km 
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POD  CSI  r— :  1  l  POD'  CSI*  rT  I""!  POD"  CSI"  far 


Radar:  CP-4,  Doppler  Radius:  6  km 
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Radar:  NOAA-C,  Doppler  Radius:  4  km 


7  dBZ  20  dBZ 


Radar:  NOAA-D,  Doppler  Radius:  2  km 
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